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Abstract 

The simplest extension of the Standard Model (SM) that generally introduces new sources 
of flavour violation and CP violation as well as right-handed (RH) currents is the addition 
of a C/(l) gauge symmetry to the SM gauge group. If the corresponding heavy gauge boson 
(Z') mediates FCNC processes in the quark sector at tree-level, these new physics (NP) 
contributions imply a pattern of deviations from SM expectations for FCNC processes that 
depends only on the couplings of Z' to fermions and on its mass. This implies stringent 
correlations between AF = 2 and AF = 1 observables which govern the landscape of the 
allowed parameter space for Z'-models. Anticipating the Flavour Precision Era (FPE) ahead 
of us we illustrate this by searching for allowed oases in this landscape assuming significantly 
smaller uncertainties in CKM and hadronic parameters than presently available. To this end 
we analyze IS.F = 2 observables in — and B^^ — B^^ systems and rare K and B decays 
including both left-handed and right-handed Z'-couplings to quarks in various combinations. 
We identify a number of correlations between various flavour observables that could test and 
distinguish these different Z' scenarios. The important role of 6 — )• si^i~ and b — )• siyi> 
transitions in these studies is emphasized. Imposing the existing flavour constraints, a rich 
pattern of deviations from the SM expectations in Bg^d and K meson systems emerges provided 
Mz' < 3TeV. While for Mz' > 5TeV Z' effects in rare Bg ^ decays are found typically below 
10% and hard to measure even in the FPE, K — )• -kvu and Kl — )• 7r^£~^£~ decays provide an 
important portal to scales beyond those explored by the LHC. We apply our formalism to NP 
scenarios with induced flavour changing neutral Z-couplings to quarks. We find that in the 
case of Bd and K decays such Z-couplings still allow for sizable departures from the SM. On 
the other hand in the Bg system, constraints on 6 — >■ si~^i~ transitions basically eliminate NP 
effects from such couplings. 
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1 Introduction 

Elementary particle physicists are eagerly waiting for clear signals of New Physics (NP) 
from the LHC While the recent discovery of a scalar particle with a mass of 125 GeV 
and the unexpectedly high direct CP violation in the charm decays could already be 
such signals, presently in both cases the SM explanations of these events are possible. 
In the first case it should be possible with increased statistic to answer the question 
whether the new particle observed at the LHC is the SM Higgs boson or another one 
belonging to a particular NP scenario. In the second case the situation is less optimistic 
in view of hadronic uncertainties but the measurements of other flavour observables in 
charm decays may tell us in due time whether the events seen by the LHCb is NP or 
not. 

After numerous proposals for the physics beyond the SM in the last 35 years it is 
really time that we know which one if any of these proposal is realized in nature. In 
particular, an exciting question is whether beyond the SM Higgs, the first new particle 
to be discovered will be a new heavy gauge boson, a new heavy fermion or a new heavy 
scalar. If this discovery is to be made directly in high energy collisions then the only 
collider in this decade that could achieve this goal is the LHC. But what if nature is not 
nice to us and the lightest new particle has a mass of 5 — 10 TeV and will just escape a 
convincing detection at the LHC. While this is fortunately only a nightmare at present 
and many new particles could still be discovered by the LHC in the coming months and 
years, we cannot presently exclude the possibility that we will have to search for new 
particles first indirectly. In such a case the high precision flavour dedicated experiments 
will be of paramount importance. However, this will require the measurements of very 
many observables and a significant reduction of hadronic uncertainties in several of them 
through improved treatment of QCD effects, in particular improved lattice calculations. 

Now over the last decades significant efforts have been made by theorists to suppress 
flavour changing neutral current (FCNC) processes so that they are absent at tree- 
level. In addition to the GIM mechanism [l] that governs the flavour physics in the 
SM, the frameworks of constrained Minimal Flavour Violation (CMFV) [2]-|4] and Min- 
imal Flavour Violation at large (MFV) [5] were very instrumental in suppressing new 
flavour and CP-violating phenomena below the present experimental bounds even in 
the presence of new particles with masses of a few hundreds GeV. Selected reviews with 
comprehensive list of references can be found in [6|[7]. 

However, if the scale of NP is shifted to 5 — 10 TeV or even higher energy scales this 
kind of suppression is less important as FCNC processes are then naturally suppressed 
by the large scales of heavy particles mediating these phenomena. In fact while loop 
diagrams, like penguin diagrams of various sorts and box diagrams dominated the 
physics of FCNC processes in the last thirty years both within the SM and several 
of its extensions, we should hope at first sight that in the case of new particles with 
masses above 10 TeV this role will be taken over by tree-level diagrams. The reason 
is simple. Internal particles with such large masses, if hidden in loop diagrams, will 
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quite generally imply very small effects in FCNC processes that will be very difficult 
to measure. On the other hand tree diagrams could still in principle provide a large 
window to these very short distance scales. 

We will demonstrate in the present paper that in the simplest extensions of the SM 
which contain just a new heavy neutral gauge boson {Z') with flavour-changing quark 
couplings, the correlations between AF — 2 and AF — 1 observables in the quark 
sector, in the absence of new heavy fermions and scalars have a significant impact on 
this optimistic expectations. In fact we find that these correlations preclude NP effects 
above 10% in rare B decay branching ratios and CP-asymmetries if Mz' > 5 TeV. Much 
larger effects are still possible in rare K decays. 

The reason is simple. A tree-level Z' contribution to AF = 2 observables depends 
quadratically on A^l^^j^{Z')/Mz', where A^-^^(Z') are flavour-violating couplings with 
i,j denoting quark flavours. For any high value of Mz', even beyond the reach of the 
LHC, it is possible to flnd couplings A^-'^(Z') which are not only consistent with the 
existing data but can even remove certain tensions found within the SM. The larger 
Mz', the larger A^"'^(Z') are allowed: A^l^^{Z') ^ aijMz> with aij sufficiently small to 
agree with AF = 2 data. Once A^l^j^{Z') are fixed in this manner, they can be used 
to predict Z' effects in AF = 1 observables. However here NP contributions to the 
amplitudes are proportional to A^-^^(Z')/M|, and with the couphngs proportional to 
Mz', Z' contributions to AF = 1 observables decrease with increasing Mz'. 

Our analysis demonstrates that for 1 TeV < Mz' < 3 TeV, still in the reach of the LHC, 

indirect Z' effects can be well tested by means of rare K and B decays. For such values 
of Mz' effects up to 50% at the level of the branching ratios and measurable effects in 
CP-asymmetries are possible for Bg^d meson systems. However for Mz' > 5 TeV, this 
begins to be very difficult even in the FPE as NP effects in rare and CP-violating B 
decays turn out to be typically below 10%. Significantly larger effects are still allowed 
in rare K decays. 

On the other hand it is evident from this discussion that flavour- violating Z-couplings, 
that arise in various extensions of the SM, could in the presence of much lower value 
of Mz provide clear NP effects in rare K and B decays even if NP generating these 
couplings is outside the reach of the LHC. In this manner flavour-violating Z couplings, 
similarly to Z' couplings in rare K decays, could turn out to be an important portal to 
short distance scales which cannot be explored by the LHC. We will demonstrate that 
this is still the case for rare and K decays but not any longer for Bg decays and 
related CP asymmetries. 

In this spirit, we will flrst ask in the present paper the following question: 

• What can be learned about NP from precise measurements of flavour observables 
to be performed in this decade if the lightest messenger of NP is a heavy Z' 
gauge boson with arbitrary couplings to quarks and arbitrary mass? In particular 
we will ask the question whether it is possible to determine all these couplings 
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entirely with the help of quark-flavour violating observables for Mz' in the reach 
of the LHC To this end we will assume that the flavour diagonal couplings of Z' 
to leptons have been determined in pure leptonic processes. 

While there is a very rich literature on FCNC processes mediated by a Z'{Z) gauge 
boson and several extensive analyses have been presented on various occasions to our 
knowledge this specific question has not been addressed so far. After having positively 
answered this question we will ask the second question: 

• What are the correlations between various flavour observables in this simple frame- 
work and how do they compare with the stringent correlations implied by the 
simplest BSM frameworks on the market, the class of models with constrained 
Minimal Flavour Violation (CMFV) [2]-[4] and the models with U{2)^ flavour sym- 
? 



metry 23-30 



The simple model analyzed here can be considered as a part of a bigger theory as 
already analyzed in numerous papers in the literature. Moreover its simplicity provides 
an analytic insight into the departure from SM expectations in flavour physics and the 
role of right-handed (RH) currents. In fact certain lessons gained from the involved 



studies of the LHT model [31 32 and RS scenario with custodial protection (RSc) |33j 
as summarized in 34 will be seen here in a much simpler setting. In particular our 
analysis in the K system, where we investigate the correlation between ek and the 
K — )■ TTuu decays, can be considered as an explicit dynamical example for the findings 



of 34 



While the violation of stringent relations of CMFV is evident in this framework due to 
new sources of flavour and CP violation, it is of interest to impose t/(2)^ symmetry on 
Z' couplings and study its phenomenological implications. Such a study is more trans- 
parent than in more complicated models in which loop diagrams with heavy fermions 
and scalars accompanied by many free parameters dominate the NP contributions to 
FCNC processes. 

The analysis of Z' flavour physics presented here can be considered as a generalization of 



our recent paper 35 in which we have analyzed in detail the pattern of flavour violating 
Z' tree-level contributions in a specific 331 model (the 331 model). The generalization 
in question is three-fold: 

• First of all we consider general structure of Z' couplings to SM quarks so that at 
the fundamental level there are no correlations between flavour violation from NP 
in K, Bd and Bs meson systems. While certain correlations between them could 
be generated once the constraints from experimental data are imposed, significant 
NP effects in and in particular in rare K decays are now possible, while this 
was not the case in the 331 model. 



^It is not possible to refer to all these papers. Selected analyses can be found in [8-21 . See also 



the review in 22 
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• Also a new important feature of NP contributions in the present paper is the pres- 
ence of flavour-violating right-handed (RH) Z' couplings to SM quarks, which has 
profound implications for correlations between AF = 2 and AF = 1 observables 
identified in scenarios with only left-handed (LH) couplings. Also correlations 
between rare decays with vv and [i^ ii~ in the final state can be modified in a 
profound manner. 

• While in the 331 model the flavour diagonal couplings of Z' to neutrinos and 
muons were fixed and smaller than the corresponding ordinary Z couplings, in 
a general case considered here they could be larger than the latter, enhancing 
thereby the branching ratios for rare leptonic and semi-leptonic decays for fixed 
quark couplings. 



As advertised above, our anatomy of Z' scenarios will lead us to the conclusion that the 
correlations between various flavour observables will test this type of NP in the FPE 
provided Mz' < 3TeV. For M^/ > 5TeV this will be very difficult, except for rare K 
decays and in the second part of our paper we will apply our formalism to the case of 
flavour-violating Z couplings. Here the effects in rare 5^ decays and in particular K 
decays can be much larger than those allowed in the case of Z' for Mz' > 1 TeV, but 
in the Bg system significant NP effects from flavour violating Z coupling are already 
ruled out by present constraints from h — )■ s[i^ y^T transitions. Similar conclusions have 



been reached in 36,37 in a more general context. 

For readers interested mainly in our results and less in the formalism presented in 
subsequent sections we have made an overview of all correlations and anticorrelations 



found by us and of the related figures in Tables [9] and 10 The comments in the last 



column of this table indicate the relevance of a given correlation or anticorrelation. 

Our paper is organized as follows. In Section |2] we describe our strategy by listing pro- 
cesses to be considered. Our analysis will only involve processes which are theoretically 
clean and have simple structure. Here we will also introduce a number of different sce- 
narios for the Z' couplings to quarks thereby reducing the number of free parameters. 
In Section [3] we will first present a compendium of formulae for master functions that 
govern FCNC processes with Z' contributions taken into account. Subsequently we 
present formulae for flavour observables in AF = 2 transitions including for the first 
time NLO QCD corrections to tree-level Z' contributions. Finally formulae for rare K 
and B decays considered by us are collected. In Section |4] we calculate Z' contributions 
to the B — )■ decay improving on the calculation of QCD corrections present in the 
Z'-literature by using the general formulae of p8] . In Section |5] we present a general 
qualitative view on NP contributions to flavour observables in four scenarios for the Z' 
couplings. In Section [6] we present our strategy for the numerical analysis and in Sec- 
tion [7] we execute our strategy for the determination of Z' couplings and discuss several 
scenarios of its couplings in question, identifying stringent correlations between various 
observables. In Section [s] we investigate what the imposition of the f/(2)^ flavour sym- 
metry on A^-'(Z') couplings would imply. In Section [o] we apply our formalism to the 
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SM Z boson for which the mass Mz and the diagonal lepton couphngs are known. A 



summary of our main results and a brief outlook for the future are given in Section 10 



2 Strategy 

Our paper is dominated by tree-level contributions to FCNC processes mediated by 
a heavy neutral gauge boson Z' . These contributions are governed by the couplings 
A^^(Z') to quarks and the corresponding Feynman rule has been shown in Fig 1 
Here j) denote quark flavours. As we will see in addition to a general form of these 
couplings it will be instructive to consider the following four scenarios for them keeping 
the pair (i, j) fixed: 

1. Left-handed Scenario (LHS) with complex A^" ^ and = 0, 

2. Right-handed Scenario (RHS) with complex A^ ^ and A^^ = 0, 

3. Left-Right symmetric Scenario (LRS) with complex A^" = Ag ^ 0, 

4. Left-Right asymmetric Scenario (ALRS) with complex A^^ = -Ag ^ 0, 

with analogous scenarios for the pair (s,(i). We will see that these simple scenarios 
will give us a profound insight into the flavour structure of models in which NP is 
dominated by left-handed currents or right-handed currents or left-handed and right- 
handed currents of the same size. In particular the last two scenarios will exhibit a very 
clear distinction between K — )■ -nvv decays and B^^d A*^/^" which are governed by 
V and A couplings, respectively. Moreover we will consider a scenario with underlying 
flavour t/(2)'^ symmetry which will imply relations between and A^^ couplings and 
interesting phenomenological consequences. 

The idea of looking at the first three NP scenarios is not new and has been in particular 
motivated by a detailed study of supersymmetric flavour models with NP dominated by 
LH currents, RH currents or equal amount of LH and RH currents 39p Moreover, it has 



been found in several studies of non-supersymmetric frameworks like LHT model 31 
or Randall-Sundrum scenario with custodial protection (RSc) [33j that models with the 
dominance of LH or RH currents exhibit quite different patterns of flavour violation. 
Our analysis will demonstrate it in a transparent manner. 

Let us then outline our strategy for the determination of Z' couplings to quarks and for 
finding correlations between flavour observables in the context of the simple scenarios 
listed above. Our strategy will only be fully effective in the second half of this decade, 
when hadronic uncertainties will be reduced and the data on various observables sig- 
nificantly improved. It involves ten steps including a number of working assumptions: 



^Similar scenarios have been considered subsequently in 36 40 
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Z' 



i^,5^p Al{Z')PL + /^tiZ')PR 



J/3 



Figure 1: Feynman rule for the coupling of a colourless neutral gauge boson Z' to quarks, 
where i, j denote different quark flavours and a, (3 the colours. Pl,r = (1 =F 75)/2- 



Step 1: 

Determination of CKM parameters by means of tree-level decays and of the necessary 
non-perturbative parameters by means of lattice calculations. This step will provide the 
results for all observables considered below within the SM as well as all non-perturbative 
parameters entering the NP contributions. As \Vub\ is presently poorly known, it will 
be interesting in the spirit of our recent papers 30 ,41 to investigate how the outcome 
of this step depends on the value of \Vub\ with direct implications for the necessary size 
of NP contributions which will be different in different observables. 

Step 2: 

We will assume that the ratios 

Mz> ' Mz' 

have been determined in pure leptonic processes. We will further assume that these 
ratios are real but could have both signs. In principle these ratios can be determined up 
to the sign from quark flavour violating processes considered below but their knowledge 
increases predictive power of our analysis. In particular the knowledge of their signs 
allows to remove certain discrete ambiguities and is crucial for the distinction between 
LHS and RHS scenarios in Bs^d ^~^fi~ decays. 

Step 3: 

Here we will consider the B^ system and the observables 

AM,, S^^, i3(5,->MV), 5^^-, (2) 



where 5'^+^- measures CP violation in Bg — )■ //"'"/i decay |42||43]. Explicit expressions 



for these observables in terms of the relevant couplings caiiDefound in Section [3j 

Concentrating in this step on the LHS scenario, NP contributions to these three ob- 
servables are fully described by 

Mz' Mz' ' Mz' ' ^ ' 

with the second ratio known from Step 2. Here S23 > and it is found to be below 
unity but it does not represent any mixing parameter as in [35]. The minus sign is 
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introduced to cancel the minus sign in Vts in the phenomenological formulae listed in 
the next section. 

Thus we have four observables to our disposal and two parameters in the quark sector to 
determine. This allows to remove certain discrete ambiguities, determine all parameters 
uniquely and predict correlations between these four observables that are characteristic 
for this scenario. 

Step 4: 

Repeating this exercise in the B'^ system we have to our disposal 

AMrf, S^Ks, i5(5rf^/iV), -^^v- (4) 

Explicit expressions for these observables in terms of the relevant couplings can be 
found in Section HI 

Now NP contributions to these three observables are fully described by 

Mz> Mz> ' Mz' ' ^ ' 

with the last one known from Step 2. Again we can determine all the couplings uniquely 
to be used in the steps below. Our notations and sign conventions are as in Step 3 with 
§13 > but no minus sign as Vtd has no such sign. 

Step 5: 

Moving to the K system we have to our disposal 

Bk, K+^n+uu, Kl^t^^v, i^L ^ 7r°£+r, fsTi ^ /iV", (6) 

where in view of hadronic uncertainties the last decay on this list will only be used 
to make sure that the existing rough bound on its branching ratio is satisfied. In the 
present paper we do not study the ratio e' /e, which is rather accurately measured but 
subject to much larger hadronic uncertainties than observables listed in ([6]). Explicit 
expressions for the observables in the K system in terms of the relevant couplings can 
be found in Section [H 

Now NP contributions to these four observables are fully described by 

Mz' Mz' ' Mz' ' Mz' ^ ' 

where we assumed the couplings to neutrinos to be left-handed and real. The ratios 
involving leptonic couplings are known already from Step 2. Consequently, we can 
determine all couplings involved by using the data on the observables in (|6]). Moreover 
we identify certain correlations characteristic for LHS scenario. Si2 > and the minus 
sign is chosen to cancel the one of Vts- 
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Step 6: As all parameters of LHS scenario has been fixed in the first five steps we are 
in the position to make predictions for the following processes 

5^x,£+r, B^Ktr, B^K*tr (8) 

B Kuu, B K*uu, B X,uu, (9) 
B X,7, B K*-f (10) 
and test whether they provide additional constraints on the couplings. 
Step 7: 

We repeat Steps 3-6 for the case of RHS. We will see that in view of the change of the 
sign of NP contribution to Bs^d — ^ fi~^fi~ and Kl — )■ fi~^fi~ decays the structure of the 
correlations between various observables will distinguish this scenario from the LHS one. 
Yet, as we will find out, by going from LHS to RHS scenario we can keep results of Steps 
3-5 unchanged by interchanging simultaneously two big oases in the parameter space 
that we encountered already in our study of the 331 model. This LH-RH invariance 
present in Steps 3-5 can be broken by the b — )■ si~^i~ and b — )■ suu transitions listed in 
([s]) and ([9]), respectively. They will allow us very clearly to distinguish the physics of 
RH currents from LH ones. As only RH couplings are present in the NP contributions 
in this scenario we can use the parametrization of these couplings as in (|3]), (|5]) and ([T]) 
keeping in mind that now RH couplings are involved. 

Step 8: 

We repeat Steps 3-6 for the case of LRS. Here the new feature is the vanishing of 
NP contributions to Bg^d — >■ ^'^ and Kl — )■ decays and rather sizable NP 

contributions to AF = 2 observables due to the presence of LR operators. As the LH 
and RH couplings are equal we can again use the parametrization of these couplings as 
in ([3]), ([5]) and ([T]) but their values will change due to different constraints from AF = 2 
transitions. Also in this step the b — )■ s£^£~ and b svv transitions will play very 
important role. 

Step 9: 

We repeat Steps 3-6 for the case of ALRS. Here the new feature is the vanishing of NP 
contributions to — )■ vr+z/P and — )■ Ti^uv decays, while Bg d including 
S^t^- CP asymmetries and again the b — )■ s(!.^(!.~ and b ^ svv transitions will exhibit 
their strength in testing the theory in a different environment: rather sizable NP con- 
tributions to AF = 2 observables due to the presence of LR operators. As the LH 
and RH couplings differ only by a sign we can again use the parametrization of these 
couplings as in ([3]), ([s]) and ([T]) but their values will change due to different constraints 
from AF = 2 transitions. 

Step 10: 

One can consider next the case of simultaneous LH and RH couplings that are unrelated 
to each other. This step is more challenging as one has more free parameters and in 
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order to reach clear cut conclusions one would need a concrete model for Z' couplings 



or a very involved numerical analysis 36 , 40 . We will therefore leave out this step from 
our paper. 

Once this analysis of Z' contributions is completed it will be straightforward to apply 
it to the case of the SM Z boson with flavour violating couplings. 

We should remark that we have left from this analysis e' je. This ratio is important 
for the tests of Z'{Z) FCNC scenarios as it is very sensitive to any NP contribution 



44-46 . However, due to significant hadronic uncertainties it is less suitable for the 
determination of the Z'{Z) FCNC couplings than the decays used by us when the 
latter will be precisely measured. On the other hand having these couplings one can 
make predictions for e' /e and study correlations. 



3 Compendium for the Contributions 
3.1 Parametrization 

First it will be useful to introduce a useful parametrization of Z' contributions by 
generalizing the master functions known from CMFV models and the LHT model to 
include in addition to left-handed currents also right-handed currents. In the case of 
the RSc model this has already been done in |33j| but our parametrization below is a 
bit different than the one in the latter paper. For our purposes it will be sufficient to 
consider the following functions: 

• For AF = 2 processes 

SiK), S{B,), S{B,), (11) 

where we will include in the definitions of these functions the contributions of 
operators with LL, RR and LR Dirac structures. 

• For decays with uu in the final state 

XL,R(i^), XL,R(5rf), Xl,r(S,). (12) 

• For decays with fift in the final state 

Ya{K), Yj,{B,), Fa(5,). (13) 



All these functions in contrast to the SM and more generally CMFV models depend on 
the meson considered and moreover are complex valued. 
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In the SM the corresponding flavour universal real valued functions are given as follows 

Axt-Ux^t+^t Sxj^logxt 

^°<"'' = 4(1 -x.r^ — ' 



Xoixt) 



Xf 



Xt + 2 3xt — 6 J 

xt-1 {xt-iy 



(15) 



^«(-') = 8 [—1 + 1)2 j • (16) 

In other CMFV models they take different values still keeping flavour universality and 
being real valued. This implies very stringent relations between various observables in 
three meson system in question which have been reviewed in ^ . It is evident that in the 
presence of Z' tree-level contributions the breakdown of flavour universality and also the 
presence of new complex phases implies the violation of these relations. Generalizing 



the three SM functions to twelve functions listed in (11)-(13), allows to describe these 
new effects in a transparent manner. In what follows we will list explicit expressions 
for these functions. Subsequently we will show how they enter the branching ratios for 
various decays. 

The derivation of the formulae listed below is so simple that we will not present it 
here. From the normalization of the corrections from Z' to the master functions in 
question and the formulae for observables given subsequently, it will be clear how these 
functions have been defined in the corresponding effective Hamiltonians. In any case, 
the compendium given below is self-contained as far as numerical analysis is concerned. 

3.2 Master Functions Including Z' Contributions 

Calculating the contributions of Z' to various decays it is straightforward to write down 
the expressions for the master functions in terms of the couplings defined in Fig. [ij 

3.2.1 AF = 2 Master Functions 

We define the relevant CKM factors 

= V*V.,, Xi'^ = VtlVt,, A? ) = VtlVts, (17) 

and introduce 

91m = ."2. = 1-78137 X lO-^CeV^^ . (18) 

V2 27r sm Ow 
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The AF = 2 master functions for M = K, Bq are given as follows 

S{M) = S^{xt) + AS(M) = \S{M)\e"'s (19) 
with AS'(M) receiving contributions from various operators so that it is useful to write 
A^(M) = [AS(M)]vLL + [AS(M)]vRR + [AS(M)]lr. (20) 



The contributing operators are defined for the K system as follows 47 , 48 







(21a) 




= (s^.PRd) (sYPRd) , 


(21b) 




= (s^^PLd) (s^PRd) , 


(21c) 




= (sPLd) (sPRd) . 


(21d) 



with analogous expressions for B^ d systems. For instance in the Bg system Q^^^ = 
(b'-y^PLs) (b'j^PLs) . Here we suppressed colour indices as they are summed up in each 
factor. For instance S'^nPid stands for Sa'^nPida and similarly for other factors. 

[A5'(M)]vLL and [AS'(M)]vrr can be obtained directly from our previous paper |35): 



[A5(5,)]vLL 



A?(Z') 



A 



(<?) 



4r 



[Asm 



VLL 



X 



4r 



(22) 



where r = 0.985 for Mz' = ITeV. [AS'(M)]vrr is then found from the formula above 
by simply replacing L by R. For the case of tree-level Z exchanges f = 1.068. 

In order to calculate the LR contributions we introduce quantities familiar from SM 
expressions for mixing amplitudes 



TiK) 



(23) 
(24) 



127r2 

where rji are QCD corrections and Bi known SM non-perturbative factors. 
Then 

TiK)[ASiK)U = [C{^^(/i^O(Q^^^(/iz',i^)) +C2^^(/iz')(Q^^(/iz',i^))] 

(25) 
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Including NLO QCD corrections 48 the Wilson coefficients of LR operators are given 
by 



2 a 
6 



47r 



—6 log 



Next 



(26) 
(27) 

(28) 



are the matrix elements of operators evaluated at the matching scale /i^' = 0{Mz> 



and are the coefficients introduced in 47 . The /i^/ dependence of C^{iJ,z') cancels 
the one of P^^jiz') so that S{K) does not depend on /i^/. 

Similarly for systems we have 

T{B,)[^S{B,)U = ^"^'^^'^^^^^'^ [^^i"^(/^z')(Q^^^(/^z',5,)) +C2^^(^z')(Q^^(/^z',i?.))] 

(29) 

where the Wilson coefficients C^{^z') are as in the K system and the matrix elements 
are given by 



(30) 



Finally, we collect in Table [ij central values of {Qtil^'Z'))- They are given in the MS-NDR 
scheme and are based on lattice calculations in 49 , 50 for — K'^ system and in jsi] 
for 5° ,, — Bds systems. For the — system we have just used the average of the 
results in 49 , 50 that are consistent with each other. As the values of the relevant Bi 
parameters in these papers have been evaluated at = 3 GeV and 4.2 GeV, respectively, 



we have used the formulae in 47 to obtain the values of the matrix elements in question 
at II z'- For simplicity we choose this scale to be Mz' but any scale of this order would 
give the same results for the physical quantities up to NNLO QCD corrections that are 
negligible at these high scales. The renormalization scheme dependence of the matrix 
elements is canceled by the one of the Wilson coefficients. 

In the case of tree- level Z-exchanges we evaluate the matrix elements at mt(mj) as the 
inclusion of NLO QCD corrections allows us to choose any scale of 0{Mz) without 
changing physical results. Then in the formulae above one should replace Mz' by Mz 
and /iz' by mt{mt). The values of hadronic matrix elements at mt{mt) in the MS-NDR 
scheme are given in Table pin 



^We thank Robert Ziegler for checking the results in this table. 
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-0.14 


0.22 


-0.11 


0.18 




-0.25 


0.34 


-0.21 


0.27 




-0.37 


0.51 


-0.30 


0.40 



Table 1: Hadronic matrix elements (Qi) in units of GeV^ at Mz' = ITeV and at 
mtimt) 



3.2.2 AF = 1 Master Functions 

We find 



Af{Z') A'j^{Z') 

X^{B^) 



AtiZ') 



AtiZ') 



[Af{Z' 



MWsM. 

^At{Z') - A^iiZ') 



Ya{B,) = r]YYo{xt) + 



[Af{Z' 



Af (ZQ - At{Z' 



qbf 



\YA{B,)\e<\ 



Here 7]x,y are QCD factors which for = mt{mt) are close to unity 52,53 

r]x = 0.994, r]Y = 1.012 . 

3.2.3 Effective Hamiltonian for b si+i' 



(31) 
(32) 

(33) 
(34) 



^\YAme''-, (35) 



(36) 



(37) 



For our discussion of constraints from b — )■ si~^i transitions we will need the corre- 
sponding effective Hamiltonian which is a generalization of the SM one: 

H efr(6 sii) = H eff(6 -> s^) - ^^V*Vt, [C^{^^)Q^i^^) + C'.iMif^)] (38) 

^ j=9,10 
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where 

Qg = (^7m^l&)(^7'^^), Qio = {s^^PLhWl^ld) (39) 

Q'g = {s^^Pnh){l^n), Q\, = {-si,PRh){lr^d) ■ (40) 

Here l-L es{b — ?■ 57) stands for the effective Hamiltonian for the b ^ s'j transition that 
involves the dipole operators. An exphcit formula for the latter Hamiltonian will be 
presented in the next section. For the Wilson coefficients we find 

1 1 7')A^^( 7') 

sin'OwC, = [vYYoixt)-Asm'ewZoixt)] - ——^ J.J > (41) 

1 1 Af(Z')A'/(Z') 



sin^^H^Cio = -VvYoixt) - -^j^ v*V ' ^^2) 

bf 7/\ A f^P-( 7') 

(43) 



sm Ow<^g - --^jV2 77^ 



•2, _ 1 1 Af(^OAT(^0 

9sM^^h' ^ts^b 



where we have defined 



(44) 



Af (Z') = A^''(Z') + Af (Z'), 
Af{Z') = Af{Z') - AfiZ'). 



(45) 



Here Zq^xi) is the SM one-loop function, analogous to Xo{xt) and Yo{xt), that represents 
gauge invariant combination of Z— and photon penguin diagrams: 

1 , 18x^ - 163x^ + 259x2 - 108a; 32x* - 38x^ - 15x^ + 18x , 
Zo(x) = -- logx + 144(^ _ 1)3 + • 

(46) 

The presence of additional coupling A'^{Z') or A^^(Z'), in addition to A'^{Z'), intro- 
duces two new parameters and allows thereby to avoid present constraints on the coef- 
ficients Cg and Cg. Therefore only the constraints on Cio and C[q from B — > K*i^i~ , 
B Ki+i- and B Xsi+i' will be relevant in the case of Z'. In the case of FCNC 
processes mediated by Z, which will be discussed in Section [9} all leptonic couplings 
are known and also the constraints on the coefficients Cg and Cg have to be taken into 
account. 

The formulae above do not include QCD renormalization group effects which influence 
only Cg and Cg. They will be taken into account in the model independent bounds on 
these coefficients in Section |9l 
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3.3 Basic Formulae for Observables 
3.3.1 AF = 2 Observables 

The Ai? = 2 mass differences are given as follows: 



G 



AM, = ^M^mBjXi'^\'FlBB,VB\S{B,)\ , (47) 

AM, = ^^M^mBMi'^\'FlBBMS{B,)\ • (48) 
The corresponding mixing induced CP-asymmetries are then given by 

S^Ks = sin(2/3 + 2^B,) , S^^ = sin(2|/3,| - 2^bJ , (49) 

where the phases (3 and /3s are defined by 

Vtd=\Vta\e-'^, = ~\Vts\e-'^\ (50) 

(3s — —1°. The new phases (pB^ are directly related to the phases of the functions 

s{B,y. 

^^B, = -9p. (51) 
Our phase conventions are as in {35] and our previous papers quoted in this work. 
For the CP-violating parameter ek and AMk we have respectively 

= . [53 (M/^)] , AMk = 25J (M,^) , (52) 

V Z[/\MK)exp 

where 

(Mi^)* = ^F^BkttikM^ [X%x, + Xh^SiK) + 2\,\tmSo{xc, Xt)] . (53) 
Here, So{xc,Xt) is a real valued one-loop box function for which explicit expression is 



given e. g. in |54 . The factors rji are QCD corrections evaluated at the NLO level 
in 



55 59]. For rji and r/s also NNLO corrections have been recently calculated 60 ,61 
Next if, = (43.51 ± 0.05)° and = 0.94 ± 0.02 p'.'es] takes into account that ^ f 
and includes long distance effects in 53(ri2) and 53(Mi2). 

In the rest of the paper, unless otherwise stated, we will assume that all four parameters 
in the CKM matrix have been determined through tree-level decays without any NP 
pollution and pollution from QCD-penguin diagrams so that their values can be used 
universally in all NP models considered by us. 
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3.3.2 Bd,s [i^^iT 



With the assumption that the CKM parameters have been determined independently 
of NP and are universal we find 



(54) 



— > 



where YA(-Bg) is given in p6[ ). 

As stressed in |42||43}[6^ when comparing the theoretical branching ratio BiBg 

with experimental data quoted by LHCb, ATLAS and CMS, a correction factor 
has to be included which takes care of AF^ effects that influence the extraction of this 
branching ratio from the data: 



B{B, fi^fi )th = r{y,) B{Bs fx^fi )cxp, r(0) = 1. 



Here 



with 



1 + A\^ys 



1 - ^ir2/. 



0.088 ±0.014. 



(55) 



(56) 



(57) 



The quantity is discussed below. 

It is a matter of choice whether the factor r(?/s) is included in the experimental branching 
ratio or in the theoretical calculation, provided r{y^ is not significantly affected by NP. 
Once it is measured, its inclusion in the experimental value, as advocated in 42 , should 



be favoured as it would have no impact on the theoretical calculations of branching 
ratios that do not depend on AP^. As in the SM and CMFV = 1 43 and the 



factor r(?/s) is universal, it is also a good idea to include this factor in experimental 
branching ratio. In this manner various CMFV relations remain intact. 

If a given model predicts significantly different from unity and the dependence 
of '"(?/s) on model parameters is large one may include this factor in the theoretical 
branching ratio: 

B[B, ;U+/i-)th 



(58) 



The branching ratios B{Bq — >■ fj,'^fj,~) are only sensitive to the absolute value of YA{Bq). 
However, as pointed out in |42||43] in the flavour precision era these decays could allow 
to get also some information on the phase of Y^lBg) and we want to investigate whether 



in the models considered this effect is significant. The authors of 43,65) provide general 
expressions for and 'S'*+^_ as functions of Wilson coefficients involved. Using these 



We follow here presentation and notations of 42 43 . 
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formulae we find in Z' models very simple formulae that reflect the fact that Z' and 
not scalar operators dominate NP contributions: 



= cos(2^^^ - li^B, 



OS 



sin(2^^= - 2(^bJ 



(59) 



Both ^fi+fj.- theoretically clean observables. 



In the formulae (59) and (61 ) we took into account new phases in the Bg — B„ mixings 



as we deal here with the mixing induced CP violation. While smaller than the phases 
of YA^Bq) their inclusion could be relevant one day. The SM phases cancel in this 
asymmetry (43 65p " 



In the SM and CMFV models 



0, r(?/,) =0.912 ±0.014 



(60) 



independently of NP parameters. 

While AFf^ is very small and yd can be set to zero, in the case of B^ — ^^fi^ one can 

651, for which we simply find 



still consider the CP asymmetry S^+^ 



S' 



sin(2^^''-2^Bj. 



(61) 



1-9 



The most recent results from LHCb read |66l 67 
B{Bs /iV") = (3.2l}|) X 10 
B{Bd-^ fx+fi-) < 9.4 X 10^1° 



B{B, fi^fi'f^ = (3.23 ± 0.27) x 10"^ (62) 

B{Bd /i+/i")^^ = (1.07 ± 0.10) X 10^^°. (63) 

We have shown also SM predictions for these observables ^8] that do not include the 
correction r{ys). If this factor is included one finds |42|[43| 

B{Bs ^ /i+/i-),s„^, = (3.5 ± 0.3) ■ 10-^ (64) 

It is this branching that should be compared in such a case with the results of LHCb 



given above. For the latest discussions of these issues see 42 , 43 , 65 , 68 



As we will see below in the Z' models considered by us 0.5 < — ^ with the 

smallest values corresponding to the largest allowed values of \S^ip\. Thus the AF^ 
effect in question varies from 5% to 9%. In view of still large experimental error we will 
approximately include this effect in the experimental branching ratio using the values 
in (60). If this is done the experimental results in (62) is reduced by 9% and we find 



B{Bs ^ 12+ fi- 



(2.9t}:t) X 10- 



(65) 



that should be compared with the SM result in (62). While the central theoretical 



value agrees very well with experiment, the large experimental error still allows for NP 



contributions. In our plots we will show the result in (65). 



^We thank Robert Knegjens and Robert Fleischer for discussion of this point. 
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3.3.3 n+n' 



Only the so-called short distance (SD) part to a dispersive contribution to — )• /i^/i^ 
can be reliably calculated. Therefore in what follows this decay will be treated only 
as an additional constraint to be sure that the rough upper bound given below is not 
violated. We have then following [69] (A = 0.226) 



B{Kl -> /i+/i")sD = 2.08 • 10-^ [P, (Yk) + A^Rt \Ya{K)\ cos P^] 



where Rt is given in (125), |V"cb| == AX and 



Y ; 



PciYi 



K) 



(66) 



(67) 



in (50). 



with Pc {Yk) = 0.113 ± 0.017 70 . Here (3 and Ps are the phases of Vtd and Vts defined 



The extraction of the short distance part from the data is subject to considerable 

S(i^i ^ /xV)sD < 2.5 ■ 10-^ (68) 



uncertainties. The most recent estimate gives 71 



to be compared with (0.8 ± 0.1) • 10"^ in the SM 70 



3.3.4 K+ TT+iyp and Kl 7r°j//> 

These are the two theoretically cleanest rare decays in quark flavour physics. Reviews 



of these two decays can be found in 72 -74 . The branching ratios for these two modes 
can be written generally as 



V 



/ReX, 



cff 



Bt{Kl ttVz/) = Kl 



V A5 



PciX) 



where 75 



V 

= (5.36 ± 0.026) • 10"^^ , = (2.31 ± 0.01) ■ 10"^° 



and 76 79 



PciX) = 0.42 ± 0.03. 
The short distance contributions are described by 

Xeff = VZVtd{XL{K) + Xn{K)) 



(69) 

(70) 

(71) 
(72) 

(73) 
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where Xl,r{K) are given in (31) and (32). 
Experimentally we have [80] 

B{K+ ^ 7r+z/z7)exp = (17.3 



+ 11.5N . 1Q- 

10. 5J 



(74) 



and the 90% C.L. upper bound [ST] 

B{Kl 7r°z/z/)exp < 2.6 ■ 10" 



(75) 

In the SM one finds |78||82) 

B{K+ vr+z/z/)sM = (8.5 ± 0.7) ■ 10"^^ , (76) 

B{Kl ^ 7rVz/)sM = (2.6 ± 0.4) ■ 10"^^ , (77) 
where the errors are dominated by CKM uncertainties. This should be compared with 



the experimental values given in (74) and (75). Clearly we have to wait for improved 
data. 



3.3.5 B {Xs, K, K*}ui> 

Following the analysis of |83], the branching ratios of the B — t- {XsiK^K*}!^!) modes 
in the presence of RH currents can be written as follows 



B{B Kvv) = B{B Kvv)sM x [1 - 27]] , 

B{B K*vv) = B{B K*vv)sM x [1 + 1.31r/] 

B{B X,uu) = B{B XsUu)sM x [1 + 0.09//] 
where we have introduced the variables 



\X^{Bs)\^ + \Xn{B, 
\VxXoixt)\^ 



V 



-Re(XL(g.)X*(E,)) 
|Xl(5.)|2 + |Xr(5,)|2 



(78) 
(79) 
(80) 

(81) 



with Xl,r{B,) defined in (33) and (34). 



Moreover the average of the K* longitudinal polarization fraction Fl also used in the 
studies of -B — 7- K*i~^i~ is a useful variable as it depends only on rj: 



(Fl) = 0.54 



(l + 2r7) 
1 + 1.3177) ■ 



(82) 



We should remark that the expressions in Eqs. (78)-(80), as well as the SM results 
in (83), refer only to the short-distance contributions to these decays. The latter are 
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obtained from the corresponding total rates subtracting the reducible long-distance 



effects pointed out in 84 



The predictions for the SM branching ratios are [83j-|85] 

E{B Kuu)sM = (3.64 ± 0.47) x 10"^ 
B{B ^ irVz/)sM = (7.2 ± 1.1) X 10"^ , 
B{B XsVu)sM = (2.7 ± 0.2) x 10"^ , 



to be compared with the experimental bounds 86 -88] 

B{B Kuu) < 1.4 X 10"^ 
B{B K*uu) < 8.0 X 10-^ 
B{B Xsvu) < 6.4 X 10"^ 



^3) 



B4) 



As e and rj can be calculated in any model by means of (81 ) the expressions given above 



can be considered as fundamental formulae for any phenomenological analysis of these 
decays and a given model can be represented by a point in the e — r] plane. Measuring 
the three branching ratios allows uniquely to determine experimentally the point (e, rj) 
and to compare with any model result. We will illustrate this for Z' scenarios. 



3.3.6 Kl 7r°£+£- 

The rare decays Kl — )■ 7T^e~^e~ and Kl — )■ 7r°yU"'"/i~ are dominated by CP-violating 
contributions. The indirect CP-violating contributions are determined by the measured 
decays Ks — >■ Tr^i'^i~ and the parameter ek in a model independent manner. It is the 
dominant contribution within the SM where one finds |89] 

B{Kl -> 7r°e+e-)sM = ^M+'.H {l.56+',il) " 10"'' , (85) 
B{Kl -> vrV+/i-)sM = 1.41^°:^^ (0.951°;^?) . lO"" , (86) 

with the values in parentheses corresponding to the destructive interference between 
directly and indirectly CP-violating contributions. The last discussion of the theoretical 
status of this interference sign can be found in 90 where the results of 91 93 are 
critically analysed. From this discussion, constructive interference seems to be favoured 
though more work is necessary. In view of significant uncertainties in the SM prediction 
we will mostly use these decays to test whether the correlations of them with Kl — )■ 
n^uu and K^ — )■ n'^uu decays can have an impact on the latter. To this end we will 
confine our analysis to the case of the constructive interference between the directly 
and indirectly CP-violating contributions. 
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The present experimental bounds 



13{Kl 7r°e+e~)cxp < 28 • 10 



-11 



94 



BiKL 7rV+/i")exp < 38 ■ 10 



-11 



195] , 
(87) 

are still by one order of magnitude larger than the SM predictions, leaving thereby 
large room for NP contributions. In fact as our numerical analysis in Sections [7] and [9] 
demonstrates, these bounds have no impact on — )■ ir'^ui? and Kl — )■ tt^uu decays 
but the present data on K'^ — )■ tt+z/z/ do not allow to reach the above bounds in the 
Z'{Z) scenarios considered. 

In the LHT model the branching ratios for both decays can be enhanced at most by a 
factor of 1.5 



31 96 . Slightly larger effects are still allowed in RSc 33 



In the LHT model, where only SM operators are present the effects of NP can be 
compactly summarised by generalisation of the real SM functions YqIxi) and ZqIxi) to 
two complex functions and Zx, respectively. As demonstrated in the context of the 
corresponding analysis within RSc 33 , also in the presence of RH currents two complex 



functions Yk and Zk are sufficient to describe jointly the SM and NP contributions. 
Consequently the LHT formulae (8.1)-(8.8) of 31 with Yk and Zk given below can 
be used to study these decays in the context of tree-level Z' and Z exchanges. The 
original papers behind these formulae can be found in (89l[9Tl|92l[97l[98] . 



Using the formulae of 33 we find 



Y, 



K 



r]YYo{xt) + 



Af{Z') 



w 



2/\^j^{Z') 
MWsM 



ts^td 



{Z') 



(88) 



(89) 



ts^td 



where is defined as in (45). These formulae with obvious changes can also be used 



for tree-level Z exchanges considered in Section [9j 

The presence of additional coupling A'^{Z') in addition to A^^(Z'), introduces as in 
B — )■ K*i~^i~ , B — )■ Kl^l~ and B — )■ Xsl^l~ two new parameters and allows thereby 
to avoid present constraints if necessary. In our analysis we will set A'^{Z') to its 
SM value. In the case of FCNC processes mediated by Z, which will be discussed in 
Section [9} all leptonic couplings are known and the predictions for Kl — )■ 7T^e~^e~ and 
Kl — )■ 7r°/i"'"/i~ are more specific. The numerical results are presented for Z' and Z 
contributions in Sections [7] and |9| respectively. 



4 B ^ Xsj Decay 
4.1 Preliminaries 



The B — )■ decay being the first loop induced B-decay determined experimentally 
has been extensively studied within the SM and its various extensions. For our cal- 
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culation of Z' contributions to the relevant Wilson coefficients very useful turned out 
to be recent study of this decay within gauged flavour models in [38]. Indeed several 
formulae of this paper could be easily adapted to our analysis. 

Let us recall that in the SM the LH structure of the W couplings to quarks requires 
the chirality ffip, necessary for b ^ s'-f transition to occur, only through the mass of 
the initial or the final state quark. Consequently the amplitude is proportional to rrih 
or rus- In contrast in models like LR models RH couplings of to quarks allow the 
chirality flip on the internal top quark line resulting in an enhancement factor rrit/mh 
of NP contribution relative to the SM one at the level of the amplitude. However, in 
the present analysis Z' contributions to -B — )■ X^j involve only SM quarks with electric 
charge —1/3 and such an enhancement is absent. Therefore we do not expect large 
corrections to 5 — )■ Xs'y from Z' exchanges, which is good as the SM agrees well with 
the data. Still it is of interest to check the size of these contributions. In doing this 
we include QCD corrections to Z' contributions at the LO using the general formulae 



of 38 , while the SM contributions are included at the NNLO level. 



Adopting the overall normalization of the SM effective Hamiltonian we have 

n,s{b ^ S7) = -^KVtb [C7,{fib)Q7, + CsG{f^b)QsG] , (90) 
where Hb = 0{mb). The dipole operators are defined as 



In writing (90) we have dropped the primed operators that are obtained from (91) by 



replacing Pr by Pl- In the SM the primed operators (RL) are suppressed by ms/mb 



relative to the ones in (90). This is not the case for RL operators but as such contribu- 
tions do not interfere with SM contributions that is dominant in any case we will neglect 
these contributions in the case of Z' as well. We have also suppressed current-current 
operators which are important for the QCD analysis. We will include these effects in 
the final formulae at the end of this section. 

The coefficients Ci{fib) are calculated from their initial values at high energy scales by 
means of renormalization group methods. We distinguish between SM quark contri- 
butions with the matching scale = 0{mt) and the Z' quark contributions with the 
matching scale iiz' = 0{Mz')- While in the LO approximation the results depend on 
the choice of the matching scale, the experience shows that taking as the matching scale 
the largest mass in the diagram appears to be a very good choice at LO. The choices 
made above follow this strategy. 

We decompose next the Wilson coefficients at the scale Hb = 0{mb) as the sum of the 
SM contribution and the Z' contributions: 

c,{^^b) = cf^{^^b) + ^cf\^^b). m 
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We recall that for the SM coefficients at fit = Olrrit) we have [xt = m^/M^^) without 
QCD corrections 

^SM/ , N _ - 2xf , -8xf - 5xf + 7xt _ ^sm/ n /Qq\ 

C7, if^t) - _ 1^ + 24(x, - 1)3 = ^""'^ ' 

In the next subsection, we summarize the results for Z' contributions to the Wilson 
coefficients of the dipole operators at the relevant matching scale iiz' = 0{Mz'). Sub- 
sequently we will present renormalization group QCD corrections to these coefficients. 
The final formula for the branching ratio for the B — )■ Xs'-j decay that includes SM and 
Z' contributions will be presented at the end of this section. 



4.2 Z' contribution without QCD Corrections 

A general analysis of neutral gauge boson contributions to -B — )■ decay has been 
presented in |38]. In addition to SM-like LL contribution from Z' we have a new LR one, 
where L (R) stands for the Pl {Pr) projector in the basic penguin diagram involving 
the s(6)-quark. 

In what follows we present the results for a contribution of a fermion / carrying electric 
charge —1/3 and having the mass ruf. This will allow us to compute the contribution 
from SM down-quarks and in the future if necessary contributions involving new heavy 
quarks. 

We ffist decompose the Wilson coefficients ACf at the ^z' scale as the sum of the 
SM-like LL contribution and a new LR one: 

Acf;(^^o = A^^cf;if,z') + A^^cf;(/iz') , 

(95) 



Adapting the formulae of 38 to our notation and denoting by / the down-quark ex- 



changed in the diagram we find 



^ (96) 

A^^Cg{^^z') = -^A^'^C^:^{^^z'). 
with Xf = m'j/Ml and summation is over the SM down-quarks. 
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For LR Wilson coefficients we find: 

^ ^'^^'"^ = -37.W>^^^ — ^» — ^"^'^ 



(97) 



A^^C^aM = -SA^''C^'{f^z'), 



with 



The summation is over SM down-quarks. 
The following properties should be noted: 

1) As opposed to the case of contributions the factor nif/mb is either C(l) or 
smaller and LR contributions are not dominant. 

2) C^q{x) is a non-vanishing monotonic function of x and takes values in the range 
[—1, —1/4] for X from to oo. 

4.3 Final Results including QCD corrections 

In order to complete the analysis of 5 — )■ we have to include QCD corrections 
which play a very important role in this decay. In the SM these corrections are known 
at the NNLO level [99j. In the LR model a complete LO analysis has been done by Cho 



and Misiak 100 and after proper modification we can use their results in our model. 



In this context the recent analyses 38 ,41 turned out to be very useful. 
We find then 

ACf;(/ife) = Kj{iiz') ACf;(/iz') + «:8(/iz') ^Ciaiiiz') + A|r^"*(/ib) . (99) 



The last contribution in ([99j) results from the mixing of new neutral current-current 
operators generated from the Z' exchange that mix with the dipole operators. The 
renormalization group analysis of this contribution is very involved but fortunately the 
LO result is known from [38]. Therefore adapting the formulae (4.16), (4.17) and (5.6) 
of this paper to our notation we find 

A^ir-\^^,) = '^^A ^'^clii^z') + Y^^Ia ^'""cii^z'), (100) 

A=L,R A=L,R 
f=u,c,t,d,s,b 



where 



(101) 
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and 



2 



Vt*sVtb 



(102) 



The diagonal couplings A^J{Z') introduce additional parameters. For our numerical 
estimate we use their SM values. 



Finally, k's are the NP magic numbers listed in Tab. [2j that is based on [38] which 
used as{Mz = 91.1876 GeV) = 0.118. They have been obtained for fit, = 2.5 GeV as 
used in the SM calculations. We add that for fj,z' = 2.5 TeV we have ^7 = 0.427 and 
K8 = 0.128. 





200 GeV 


1 TeV 


5 TeV 


10 TeV 


Mz 


Ks 


0.524 
0.118 


0.457 
0.125 


0.408 
0.129 


0.390 
0.130 


0.566 
0.111 


u,c 
'^LL 

s,b 
'^LL 


0.039 
-0.002 
-0.040 
0.087 
0.128 


0.057 
-0.003 
-0.057 
0.090 
0.147 


0.076 
-0.002 
-0.072 
0.090 
0.163 


0.084 
-0.001 
-0.079 
0.090 
0.168 


0.030 

-0.032 
0.084 
0.116 


u,c 

^LR 

f^LR 

s,b 
^LR 
r^d 


0.085 
0.004 
-0.015 
-0.078 
0.473 


0.128 
0.012 
-0.025 
-0.092 
0.665 


0.173 
0.023 
-0.036 
-0.106 
0.865 


0.193 
0.028 
-0.041 
-0.111 
0.953 


0.065 

-0.011 
-0.070 
0.383 



Table 2: The NP magic numbers relevant for QCD calculations [38]. For completeness, 
in the last column the case of a flavour-violating Z is included. 



Using these formulae we find for Mz' = 1 TeV 

ACf;M = O(10-^). (103) 
While due to the presence of RH couplings, this contribution is by one order of magni- 



tude larger than found in the 331 model 35 , it is still negligible when compared with 



the SM value of —0.353. Therefore we will not consider B — )■ decay further. 



5 General Structure of New Physics Contributions 
5.1 Preliminaries 

We have seen in Section |2] that the small number of free parameters in each of LHS, 
RHS, LRS and ALRS scenarios allows to expect definite correlations between flavour 
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observables in each step of the strategy outhned there. These expectations will be 
confirmed through the numerical analysis below but it is instructive to develop first a 
qualitative general view on NP contributions in different scenarios before entering the 
details. 

First, it should be realized that the confrontation of correlations in question with future 
precise data will not only depend on the size of theoretical, parametric and experimental 
uncertainties, but also in an important manner on the size of allowed deviations from SM 
expectations. The latter deviations are presently constrained dominantly by AF = 2 
observables and B Xg'j decay. But as already demonstrated in |35||36||40] after 
the new data from the LHCb, ATLAS and CMS also the decays Bg fi — )■ and 
b — si~^£~ begin to play important roles in this context. We will see their impact on 
our analysis as well. 

Now, in general NP scenarios in which there are many free parameters, it is possible with 
the help of some amount of fine-tuning to satisfy constraints from AF = 2 processes 
without a large impact on the size of NP contributions to AF = 1 processes. However, 
in the Z' scenarios considered here, in which NP in both AF = 2 and AF = 1 processes 
is governed by tree-diagrams, the situation is different. Indeed, due to the property of 
factorization of decay amplitudes into vertices and the propagator at the tree-level, 
the same quark flavour violating couplings and the same mass Mz' enter AF = 2 and 
AF = 1 processes undisturbed by the presence of fermions entering the usual box and 
penguin diagrams. Let us exhibit these correlations in explicit terms. 



5.2 AF = 1 vs. AF = 2 Correlations 

In order to obtain transparent expressions we introduce 



,VLL 



„VRR 



9sM 



(104) 



which is the same for K and Bg systems, f 1 is defined in (22). In the case of LR 



contributions to AS* let us rewrite (25) and (29) as follows 



[AS{K)U 



„LR/ 



(K) Af{Z')A^i{Z' 



[A 



Wl2 



(105) 



[A5(i?,)]LR 



^LR/ 



(5,) A'l{Z')A'^{Z' 



bq. 



Ml 



(106) 



where the quantities r can be found by comparing these expressions with (25) and 



(29), respectively. They depend on low energy parameters, in particular on the meson 



system and logarithmically on Mz^ The latter dependence can be neglected for all 
practical purposes as long as Mz' is in the range of a few TeV. 
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We can then derive the following relations between shifts in the basic functions in 
AF = 1 and AF = 2 processes which are independent of any parameters like Sij but 
depend sensitively on Mz' |5 In particular they do not depend explicitly on whether 
SI or S2 scenarios for \Vub\ are considered. This dependence is hidden in the allowed 
shifts in AS{K) and AS{Bd) both in magnitudes and phases. We have then 

LHS Scenario 

AXUB^ = = 0.597 (108) 



and 



AY^{K) = ^X,{K)^^^, AY^{B,) = A^l(5,)^|^. (109) 



RHS Scenario 



AXniK) = AXl{K) = (HO) 



AXniB,) = AX,{B,) = ~^YAiB,)^^^. (Ill) 



LRS Scenario 

AXl{K) AXr{K) Af{Z' 



-ASiK) ./AS{K) Mz.gl^^-r^'^'^-r^''{K) 
AXl{B,) AXn{B,) Af{Z') 



0.048, (112) 
(113) 



^-AS{B,y y/-AS{B,y Mz^gl^^-r^'^'^-r^'^iB,) 

with ^^'^^'^ = 0.204 and ^^'^^'^ = 0.212. 
^-AS{B,y V-A5(5,)* 

There are no NP contributions to functions in this scenario. 

ALRS Scenario 

AYa(K) Af(Z') , , 

' - 2 ' = = 0.094, 114 



./AS{K) Mz^gl^^r-^^^-r^^iK) 



^Similar relations have been derived in |35) in the context of 331 model but they involved only LHS 
scenario. 

^The numerical values on the r.h.s of these equations correspond to Mz' = 1 TeV. 
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^YAiB,) K'{Z') ^^^^^ 



with = 0.337 and ^dSL = 0.346. 

There are no NP contributions to Xl^r functions in this scenario. 
General Scenario 

Finally we give for completeness general formulae for the correlations in question that 
do not assume any particular relation between LH and RH couplings. To this end we 
write 

A'^ = aijA'l, aji = a*j, (116) 
where are complex numbers. 
We find then in the K system 

AXl{K) V2Af{Z') 1 



y/AS{K) Mz'qIu y/r^^^il + al) + 2a,ar^^{Ky ^^^^^ 
AXR{K)^asdAXL{K), (118) 



AY^iK) _ V2Af{Z') 1 



Similarly in the Bs,d systems we have 

AXl{B,) v^Af(Z') 1 



\/^S{B^y Mz'9iu ^/rVLL(l + a^^) + 2a,ferLR(5, 



(119) 



(120) 



AXniBg) = agbAXL{Bg) (121 
AY^iBg) _ V2A>f{Z') 1 - 



(122) 



Evidently these general relations involve more free parameters than in the scenarios 
considered in our paper but they could turn out to be useful in concrete Z' models and 
models with tree-level FCNC's mediated by Z boson. 
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5.3 Implications 

Inspecting these formulae we observe that if the SM prediction for Ek is very close to 
its experimental value AS{K) cannot be large and consequently at first sight the shifts 
AXl,r{K) and AYa{K) cannot be large implying suppressed NP contributions to rare 
K decays unless Z' couplings to neutrinos and charged leptons in the final state are 
enhanced. The details depend on the value of Mz'- However, as we will find below, the 
present theoretical and parametric uncertainties in Ek and AM^ still allow for large 
effects in rare K decays both in SI and S2 scenarios. 

Similarly in the and Bg systems if the SM predictions for AM^^d, S^Ks ^md S^^f) 
are very close to the data, it is unlikely that large NP contributions to rare Bd and Bg 
decays, in particular the asymmetries '+ , _ , will be found, unless again Z' couplings 
to neutrinos and charged leptons in the final state are enhanced. Here the situation 
concerning theoretical and parametric uncertainties is better than in the K system and 
the presence of several additional constraints from 6 — )■ s transitions allows to reach in 
the Bg system clear cut conclusions. 

In this context it is fortunate that within the SM there appears to be a tension between 
the values of Ek and S^Ks so that some action from NP is required. Moreover, parallel 
to this tension, the values of \Vub\ extracted from inclusive and exclusive decays differ 



significantly from each other. For a recent review see 101 



If one does not average the inclusive and exclusive values of \Vub\ and takes into account 
the tensions mentioned above, one is lead naturally to two scenarios for NP: 

• Exclusive (small) \Vub\ Scenario 1: \ek\ is smaller than its experimental de- 
termination, while S^Ks is rather close to the central experimental value. 

• Inclusive (large) \Vub\ Scenario 2: \ek\ is consistent with its experimental 
determination, while S^Ks is significantly higher than its experimental value. 

Thus dependently which scenario is considered we need either constructive NP con- 
tributions to \ek\ (Scenario 1) or destructive NP contributions to S^Ks (Scenario 2). 
However this NP should not spoil the agreement with the data for S^Ks (Scenario 1) 
and for \ek\ (Scenario 2). 

While introducing these two scenarios, one should emphasize the following difference 
between them. In Scenario 1, the central value of \ek\ is visibly smaller than the very 
precise data but the still significant parametric uncertainty due to \Vcb\^ dependence in 
l^j^l and a large uncertainty in the charm contribution found at the NNLO level in 61 
does not make this problem as pronounced as this is the case of Scenario 2, where large 
\Vub\ implies definitely a value of S^Ks ^h^^ is by 3cr above the data. 

Our previous discussion allows to expect larger NP effects in rare B^, decays in scenario 
S2 than in SI. This will be indeed confirmed by our numerical analysis. In the K 
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system one would expect larger NP effects in scenario SI than S2 but the present 
uncertainties in and AM^ do not allow to see this clearly. The system is not 
affected by the choice of these scenarios and in fact our results in SI and S2 are basically 
indistinguishable from each other as long as there is no correlation with the Bet system. 
However, we will demonstrate that the imposition of f/(2)^ symmetry on Z' couplings 
will introduce such correlation with interesting implications for the Bg system. 

We do not include B ^ r^Ur in this discussion as NP related to this decay has nothing 
to do with Z'. Moreover, the disagreement of the data with the SM in this case softened 
significantly with the new data from Belle Collaboration |102| . The new world average 



IS 



provided by the UTfit collaboration of B{B+ r+i/)exp = (0.99 ± 0.25) x 10"^ (TO^ 
in perfect agreement with the SM in scenario S2 and only by 1.5a above the SM value 
in scenario SI. 

Evidently \Vub\ could be some average between the inclusive and exclusive values, in 
which significant NP effects will be in principle allowed simultaneously in K and Bd 
decays. This discussion shows how important is the determination of the value of \Vub\- 

As already remarked above, the case of Bg mesons is different as the B^ — B^ system 
is not involved in the tensions discussed above. Here the visible deviation of the AM^ 
in the SM from the data and the asymmetry S^^, still being not accurately measured, 
govern the possible size of NP contributions in rare decays. 

With this general picture in mind we can now proceed to numerical analysis. 



6 Strategy for Numerical Analysis 



6.1 Preliminaries 



Similarly to our analysis in 35 it is not the goal of the next section to present a 
full-fledged numerical analysis of all correlations including present theoretical, para- 
metric and experimental uncertainties as this would only wash out the effects we want 
to emphasize. Yet, these uncertainties will be significantly reduced in the coming 
years 104, 105 and it is of interest to ask how the Z' scenarios considered here would 
face precision flavour data and the reduction of hadronic and CKM uncertainties. In this 
respect as emphasized above correlations between various observables are very impor- 
tant and we would like to exhibit these correlations by assuming reduced uncertainties 
in question. This strategy will also be used for the case of flavour violating Z-couplings. 

Therefore, in our numerical analysis we will choose as nominal values for three out of 
four CKM parameters: 

IK.I = 0.2252, iKfol = 0.0406, 7 = 68°, (123) 

and instead of taking into account their uncertainties directly, we will take them effec- 
tively at a reduced level by increasing the experimental uncertainties in AMg ^ and Ek- 
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Gf = 1.16637(1) X 10- 
Mw = 80.385(15) GeV 
sin^ Ow = 0.23116(13) 
a{Mz) = 1/127.9 
as{Mz) = 0.1184(7) 



GeV" 



m„(2 GeV) = (2.1 ± 0.1) MeV 
md{2 GeV) = (4.73 ± 0.12) MeV 

m,(2 GeV) = (93.4 ± 1.1) MeV 



mc(mc) = (1.279 ± 0.013) GeV 

m,{m,) = A.19t''om GeV 
mt{mt) = 163(1) GeV [l07 
Mt = 172.9 ±0.6 ±0.9 GeV 
rriK = 497.614(24) MeV 
Fk = 156.1(11) MeV 
Bk = 0.767(10) 

= 0.94(2) 
r/i = 1.87(76) 

= 0.5765(65) 
r/3 = 0.496(47) 

AMk = 0.5292(9) x lO-^---^ 



106 



106 



106 



106 



106 



107 



107 



107 



108 



106 



109 



106 
106 



107 



62 



07 



ps 



IekI = 2.228(11) X 10 



-3 



63 

6T 
58 
Bo 



106 



106 



ruB^ = 5279.5(3) MeV 
me, = 5366.3(6) MeV 

= (190.6 ± 4.6) MeV 
Fb, = (227.7 ±6.2) MeV 
Bb, = 1.26(11) 
Bb, = 1.33(6) 
BbJBb, = 1.05(7) 



Bb^ = 226(13) MeV 



Fb^^Bb^= 279(13) MeV 

^ = 1.237(32) 

riB = 0.55(1) 

AMd = 0.507(4) ps-i 

AM, = 17.73(5) ps-i 

S^Ks = 0.679(20) 

S^^ = 0.0002 ± 0.087 



106 



106 



107 



107 



107 



107 



107 



107 



107 



58 



110 



B{B+ 

rB± = 
\Vus\ -- 

\Vcb\ = 
I T Aincl. 
I ^ub 
I T Aexcl. 



(1.64 ±0.34) X 10" 



(1641 ±8) X 10-3 ps 

= 0.2252(9) 

: (40.6 ± 1.3) X 10-3 
= (4.27 ±0.38) X 10-3 
= (3.12 ±0.26) X 10-3 



07 



59 



06 



111 



106 



112 



106 



106 



106 



106 



106 



107 



Table 3: Values of the experimental and theoretical quantities used as input parameters. 



Here the values for \Vus\ and \Vcb\ have been measured in tree level decays. The value 
for 7 is consistent with CKM fits and as the ratio AM^/AMg in the SM agrees well 
with the data, this choice is a legitimate one. Other inputs are collected in Table |3] 
For \Vub\ we will use as two values 

iKfel = 3.1 ■ 10-3 \Vub\ = 4.0 ■ 10-3 (124) 

corresponding to central values of exclusive and inclusive determinations of this CKM 
element and representing thereby SI and S2 scenarios, respectively. 



Having fixed the three parameters of the CKM matrix to the values in ( 123 ), for a given 
\Vub\ the "true" values of the angle /3 and of the element \Vtd\ are obtained from the 
unitarity of the CKM matrix: 



\Vtd\ = \Vusmb\Ru Rt = \/l + Rl-2RbCos^ , cot/3= ^ Rbcos^ ^ 

V Kb sm 7 

where 
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Scenario 1: 


Scenario 2: 


bxperiment 




1.72(22) ■ 10-3 


2.15(32) ■ 10"3 


2.228(11) X 10-3 


(sin 2/3) true 


0.623(25) 


0.770(23) 


0.679(20) 


AM, [ps-i] 


19.0(21) 


19.0(21) 


17.73(5) 


AM, [ps-i] 


0.56(6) 


0.56(6) 


0.507(4) 




0.62(14) ■ 10-^ 


1.02(20) ■ 10-^ 


0.99(25) X 10-4 



Table 4: SM prediction for various observables for \Vub\ = 3.1-10 ^ and \Vub\ = 4.0-10 
and 7 = 68° compared to experiment. 



In Table |4| we summarize for completeness the SM results for Iek], AM,^,, (sin2/3) 



true 



and B{B^ — > r+i/^), obtained from (125), setting 7 = 68° and choosing the two values 
for \Vub\ in ( |124 ). We observe that for both choices of \Vub\ the data show significant 
deviations from the SM predictions but the character of the NP which could cure these 
tensions depends on the choice of \Vub\ as already discussed in detail in [t] and in the 
previous section. 

What is striking in this table is that the predicted central values of AM, and AM,, 
although slightly above the data, are both in good agreement with the latter when 
hadronic uncertainties are taken into account. In particular the central value of the 
ratio AMg/ AMd is very close to the data: 

fAMs\ , , 

(aI/) =34.5 ±3.0 exp: 35.0 ± 0.3 . (127) 



SM 



These results depend on the lattice input and in the case of AM, on the value of 7. 
Therefore to get a better insight both lattice input and the tree level determination of 
7 have to improve. 

In 35 we have analyzed a particular 331 model, the so-called 331 model. Because 
of suppressed contributions to Ek-, this model favoured the inclusive value of |Vu6|. 
Moreover only left-handed couplings of Z' to quarks were present. As already described 



in Section[2jthe present analysis can be considered as the generalization of 35 to include 
also exclusive values of and the right-handed couplings of Z' to quarks. Thus with 
two scenarios for \yub\ and four scenarios LHS, RHS, LRS, ALRS for flavour violating 
couplings of Z' to quarks we are led to eight scenarios of Z'-physics to be denoted by 

LHSl, LHS2, RHSl, RHS2, LRSl, LRS2, ALRSl, ALRS2 (128) 

with SI and S2 indicating the scenarios. 

We should emphasize that in each case we have only two free parameters describing 
the Z'-quark couplings in each meson system except for the universal Mz' ■ Therefore, 
as in the case of the 331 model it is possible to determine these couplings from flavour 
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observables (see Section [2]) provided flavour conserving Z' couplings to neutrinos and 
muons are known. This was the case of the 331 model. Here these couplings are 
not fixed by the theory and have to be determined in purely leptonic processes. In 
principle one could also get some insight about them from semi-leptonic meson decays 
but determining them in purely leptonic processes increases the predictive power of the 
theory. 

Following Step 2 of our general strategy of Section [2| in what follows we will assume 
that A'^{Z') and A'£'{Z') have been determined in purely leptonic processes. For 
definiteness we set the lepton couplings at the following values 

Af (Z') = 0.5, Af{Z') = 0.5, (129) 

In 



to be compared with Af{Z') = 0.14 and Af{Z') = -0.26 in the 331 model 
the SM both couplings of Z are equal to 0.372. 



35 



The specification of signs in (129) is crucial for the identification of various enhance- 
ments and suppressions with respect to SM branching ratios and CP asymmetries and 
is at the basis of our search for successful oases in the space of parameters. If these signs 
will be identified in the future to be different from the ones assumed here, it will be 
straightforward to find out by inspecting our results how the landscape of oases changes 
for each of the four possibilities for the signs of leptonic couplings. 



6.2 Dependence on Mz' 



The correlations between AF = 1 and AF = 2 derived in subsection |5.2| imply that 
when free NP parameters have been bounded by AF = 2 constraints, the modifications 
of AXi and AYi are inversely proportional to Mz'- This means that in the case of 
NP contributions significantly smaller than the SM contributions, the modifications of 
rare decay branching ratios due to NP will be governed by the interference of SM and 
NP contributions and consequently will also be inversely proportional to Mz'- This 
is the case of all observables in Bg and systems, but not in K system where NP 
contributions could be much larger than the SM contribution for sufficiently low values 
of Mz'- In the latter case the NP modifications of branching ratios will decrease faster 
than 1/Mz' (1/M|, in the limit of full NP dominance) until NP contributions are 
sufficiently small so that the 1/Mz' dependence is again valid. 

Concerning the direct lower bound on Mz' from collider experiments, the most stringent 



bounds are provided by CMS experiment 113 . The precise value depends on the model 



considered. While for the so-called sequential Z' the lower bound for Mz' is in the 
ballpark of 2.5 TeV, in other models values as low as 1 TeV are still possible. In order 
to cover large set of models, we will choose as our nominal value Mz' = 1 TeV. With the 



help of the formulae in subsection 5.2 it should be possible to estimate approximately, 
how our results would change for 1 TeV < Mz' < 3 TeV. For much larger values of 
Mz', considered mainly in K physics, explicit results will be provided. 
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6.3 Simplified Analysis 

As in j35,| we will perform a simplified analysis of Ek^ AM^ ,,5 S^Ks ^^^d S^^ in order 
to identify oases in the space of new parameters (see Section |2| for which these five 
observables are consistent with experiment. To this end we set all other input param- 
eters at their central values but in order to take partially hadronic and experimental 



uncertainties into account we require the theory in each of the eight scenarios in ( 128 ) to 
reproduce the data for sk within ±10%, AM^^^ within ±5% and the data on S^Ks ^'^^ 
S^^ within experimental 2a. We choose larger uncertainty for Ek than AM^ ^ because 
of its strong |V"cfe|^ dependence. For AMx we will only require the agreement within 
±25% because of potential long distance uncertainties. 

Specifically, our search is governed by the following allowed ranges: 

16.9/ps < AM, < 18.7/ps, -0.18 < S^^ < 0.18, (130) 
0.48/ps < AMd < 0.53/ps, 0.64 < S^Ks < 0-72. (131) 

0.75 < , . . < 1.25, 2.0 X 10"^ < {skI < 2.5 x 10^1 (132) 
(AMx)sM 



The search for these oases in each of the scenarios in ( 128[ ) is simplified by the fact that 



for fixed Mz' each of the pairs (AM,, S^^), (AM^, S^Ks) (AM^, \£k\) depend only 
on two variables. The fact that in the K system we have only one powerful constraint at 
present is rather unfortunate. The situation will improve by much when the branching 
ratios for — )■ tt'^uu and Kl — )■ n^uu will be measured. 

In what follows we will first for each scenario identify the allowed oases. As in the case 



of the 331 model there will be several oases allowed by the constraints in (130)-(132) 
and we will have to invoke other observables, which are experimentally only weakly 
bounded at present in order to find the optimal oasis in each case. Yet, our plots will 
show that once these observables will be measured precisely one day not only unique 
oasis in the parameter space will be identified but the specific correlations in this oasis 
will provide a powerful test of the Z' scenarios. 

As in |35j, inspecting the expressions for various observables in different oases, we have 
identified the fastest route to the optimal oasis in each scenario. We will describe 
this route in each case below. We will also see how the correlations between various 
observables can give additional tests once the analysis is confined to a particular oasis. 

It turns out that considerable progress in the search for the optimal oasis in each scenario 
can be made by identifying some special observables for whom the sign of departure 
from SM expectations is sufficient to identify this oasis uniquely. For instance in the 
case of LHS scenario, in the and Bg meson systems these special observables turn 
out to be 

i3(5,^/iV), (133) 
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respectively. Already the sign of shifts of them with respect to the SM values allows 
to make significant progress towards the identification of the optimal oasis in each 



scenario considered. However, in contrast to LHS scenario considered in 35 , in the 
presence of RH currents the two observables will not be sufficient to identify optimal 
oasis. As already advertised, in the case of the system, the rescue will come from 
B —7- K*^^^^ , B —7- Kfi^^^ and b — )■ sui? transitions. 



7 An Excursion through Z' Scenarios 

7.1 The LHSl and LHS2 Scenarios 
7.1.1 The Bs Meson System 

We begin the search for the oases with the Bg system as here the choice of \Vub\ is 
immaterial and the results for LHSl and LHS2 scenarios are almost identical. Basically 
only the asymmetry within the SM and \Vts\ are slightly modified because of the 
unitarity of the CKM matrix. But this changes S^^ in the SM from 0.032 to 0.042 and 
can be neglected. 

The result of this search for Mz> = 1 TeV is shown in Fig. [2| where we show the allowed 
ranges for (523,^23)- The red regions correspond to the allowed ranges for AMg, while 
the blue ones to the corresponding ranges for S'^^. The overlap between red and blue 
regions identifies the oases we were looking for. We observe that the requirement of 
suppression of AMg implies §23 7^ 0. 

From these plots we extract several oases that are collected in Table [5] We denote by 
Ai{Sl) and Ai{S2) the oases found for the two values of \Vub\ but as in the -B^-case 
there is no change in these oases when moving from SI to S2 we will show the results 
only for LHSl scenario. We observe the following pattern: 

• For each oasis with a given ^23 there is another oasis with ^23 shifted by 180° 
but the range for S23 is unchanged. This discrete ambiguity results from the fact 
that AMg and 5*^0 are governed by 2^23. However, as already seen in Table [5] 
and discussed below this ambiguity can be resolved by other observables. In this 
context we just investigate whether in a given oasis various branching ratios are 
enhanced or suppressed with respect to the SM or CP asymmetries modified. In 
the case of 5'^+^-, that vanishes within the SM, we just look at its sign. In the 
last two columns of Table [5] we consider 

ABf^~ = AB{B, fi^fi-), ABf = AB{B ^ X,uu) . (134) 

• The oases with i = 2,4 are very small and imply very concrete predictions for 
various observables. In fact as we will soon see they are already ruled out by the 
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present data on E{Bs — ^ /^^Ai )• They correspond roughly to NP contribution to 
Mfg twice as large as the SM one but carrying opposite sign. 

• The increase of Mz' by a given factor allows to increase S23 by the same factor. 
This structure is evident from the formulae for IS.S[B^. However, the inspection 
of the formulae for AF = 1 transitions shows that this change will have impact 
on rare decays, making the NP effects in them with increased M^/ smaller. This 
is evident from the correlations derived in Section |5] and has been emphasized at 
the beginning of our paper. 



We will next confine our numerical analysis to these oases, investigating whether some 
of them can be excluded by other constraints and studying correlations between various 



observables. To this end we set the lepton couplings as given in (129). 



As a final comment, we observe that the oases reported in Table |5]and all other tables 
for other scenarios below actually describe squares in the spaces (§235 (^23)5 (§135 '^13) 
and (§12, 5i2), while the corresponding regions in Fig. [2] and analogous figures for other 
scenarios have more complicated shapes. Indeed, in our numerical analysis of the various 
observables we have varied the parameters in the true oases, requiring that constraints 
are satisfied. 

In Fig. 3] (left) we show 5*^+^- vs S^^. The requirement of suppression of AMg requires 
S^+^- to be non-zero. A positive value of 'S'^+^_ chooses scenario Ai, while a negative 
one scenario ^3. Note that in both scenarios the sign of 5*^,^ is not fixed yet but it will 
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■§23 


(^23 


OS 










0.0016 - 0.0061 


49° - 129° 


+ 


± 








0.0176-0.0181 


87° - 92° 












0.0016 - 0.0061 
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Table 5: Oases in the space (523,^23) for Mz' = 1 TeV in LHSl. The sign of S^+^- 
chooses the oasis uniquely. The same applies to the pair S^^ and AB{Bs — )■ yu^/i^) as 
discussed in the text. Basically the same results are obtained in LHS2. 



LHS1 LHS1 




Figure 3: S^^+^- versus S^^ (le-ft) and S^^ versus B{Bs — > fi^fi ) (right) for Mz' = 
1 TeV in LHSl. Ai: blue, A3: purple, A^: red, A4: gray. Gray region: exp la range 
B{Bs = (2.9^};i) ■ 10^^ Red point: SM central value. 



be fixed by invoking B{Bs — ?■ fJ'^f^^) below. We note that in big oases for Mz> = 1 TeV, 
|S'^+^_| can reach values as high as 0.9 when \S^^\ ~ 0.2. Smaller values are found for 
larger Mz'- We also observe that the small oases represented by gray and red areas are 
indeed very small and imply |S'*+^_| ~ 1. 

The fact that 5"^+^- is very powerful in identifying the optimal oasis can be understood 
as follows. 5"^+^- is governed by the phase of the function Y{Bs) that originates in the 
Z' contribution. It can distinguish between Ai and ^3 oasis because the new phase ^23 
in these two oases differs by 180° and consequently sin 523 relevant for this asymmetry 
differs by sign in these two oases. Calculating the imaginary part of Y{Bs) in (36) 
and taking into account that it is Af* and not A^^** that enters Y{Bs) one can convince 
oneself about the definite sign of Sf,+ ,,. in Ai and ^3 oases as stated above. 

The reason why B{Bs — fi^fJ'") cannot be presently powerful in the search for oases is 
the significant experimental error on S^^ with which this branching ratio is correlated. 
However, inspecting this correlation in a given oasis constitutes an important test of the 
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model. We show this in Fig. [3j (right While in the oasis A\ S^^ increases (decreases) 
uniquely with decreasing (increasing) B{Bs — ?■ fi^fi~), in the oasis ^3, the increase of 
S^^ implies uniquely an increase of B{Bs — )■ fi^fi^). Therefore, while B{Bs — )■ /i"*"/^") 
alone cannot uniquely determine the optimal oasis, it can do in collaboration with S^^p. 
Finding both these observables above or below their SM expectations, would select the 
oasis A3, while finding one of them enhanced and the other suppressed (opposite sign 
in the case of S^,/,) would select Ai as the optimal oasis. We indicate this pattern in 
Table [sj In fact in the coming years it will be S^^ and B{Bs — )■ /i'''/i~) which will be 
leading this search as 3^^+ _ is much harder to measure. 

If the favoured oasis will be found to differ from the one found by means of _ one 
day the model in question will be in trouble. Indeed, let us assume that B{Bs fi ) 
will be found below its SM value. Then the measurement of S'^^ will uniquely tell us 
whether Ai or A3 is the optimal scenario and consequently as seen in Fig. |3] (left) and 



Table |5| we will be able to predict the sign of S^+^. . Moreover, in the case of 5*^ 



sufficiently different from zero, we will be able to determine not only the sign but also 
the magnitude of Sf,+ ,-. 

Probably the most important message from Fig. |3] (right) is the following one. If 
NP is dominated by Z' in the LHSl scenario, then departure of S^rp from the SM 
implies automatically the departure of B{Bs — ?■ /U^/U^) from its SM value and vice 
versa. Moreover, for Mz' = 1 TeV and \S^^\ ~ 0.2, B{Bs fJ'^fJ' ) can deviate from 
the SM value by ±60%. We also note that the small oases are inconsistent with the 
LHCb data on B{Bs — )■ /i^/x~) and are already ruled out. Consequently \S^^+^-\ ~ 1 
is also ruled out. Therefore we will omitt the results for small oases in the subsequent 
plots for Bg meson system. 

In Fig. [4] (left) we plot vs S^^. We observe that for Mz' = 1 TeV and S^tf> 

significantly different from zero, '^^^ differ significantly from unity. With 
low as 0.6 the effect of AP^ on B{Bs — t- /i^/U^) becomes smaller. 

In Fig. [4] (right) we show B{B Xsvv) vs B{Bs — t- fi^fj,^). This correlation is valid in 



any oasis due to the assumed equal sign of the leptonic couplings in (129). However, as 
seen in the plot the size of NP contribution may depend on the oasis considered. We 
note that NP effects of 50% are still possible and suppression of B{Bs /W'^/W") below 
the SM value will also imply the suppression of B{B — )■ XgUu). Yet, one should note 
that if the future data will disagree with this pattern, the rescue could come from the 
flip of the signs in uu or fi~^fj^~ couplings provided this is allowed by leptonic decays of 
Z'. 



In Fig. [5jwe show B{B — )■ XgUu) vs S^^f, which could turn out to be informative when 
S^^ will be measured precisely one day. 



« The central values for B{Bd ^i+n~f^ = 1.0 x 10"!° and B{Bs fj.+fj.-f^ = 3.1 x lO'^ 
shown in the plots correspond to fixed CKM parameters chosen by us and differ from the ones listed 



in ( 62 1 and ( 63 ) but are fully consistent with them. 
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, . ,. 3 1-6| 
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S** Bie,^;;,^:) 110-^1 

Figure 4: versus S^^ (left) and B{B — )■ Xguu) versus B{Bs l^^l^ ) (right) for 

Mz' = 1 TeV in LHSl. Ai: blue, A^: purple. Gray region: exp la range B{Bs — )■ 
= {2.9t\'f) ■ 10^^. Red point: SM central value. 
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Figure 5: B{B Xsw) versus B{Bs fi+fi') for Mz> = 1 TeV m LHSl. Ai: blue, 
As.- purple, A2: red, A^: gray. Red point: SM central value. 



7.1.2 The Bd Meson System 

We begin by searching for the allowed oases in this case. The result is shown in Fig. [6] 
and Table [6] The general structure of the discrete ambiguities is as in Table |5] but now 
as expected the selected oases in SI and S2 differ significantly from each other. 

Let us first concentrate on S2 scenario that corresponds to the one already analyzed 
in 



35 . In the right panel of Fig. [7|we show S^Ks vs B{Bd — )■ The requirement 

on S^Ks and AMa forces B{Bd — )■ fi^fi~) to differ from the SM value but the sign 
of this departure depends on the oasis considered. Here distinction is made between 
Bi and B3 for which B{Bd — )■ /^'''/i") is suppressed and enhanced with respect to 
the SM, respectively. These enhancements and suppressions amount up to ±50% for 
Mz' = 1 TeV. They increase with decreasing S^Ks- 

Note that because of the correlation between B{Bd — t- fJ'~^fJ'~) and S^pxs the fact that 
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Figure 6: Ranges for AM^ (red region) and S^Ks (blue region) for Mz' = 1 TeV 



LHSl (left) and LHS2 (right) satisfying the bounds in Eq. (131) 
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Figure 7: S^Ks versus B{Bd — )■ ) for Mz' = 1 TeV in LHSl (left) and LHS2 

(right). Bi: yellow, B^: green, B2: magenta, B^: brown. Red point: SM central value. 



the latter is aheady well determined, the range of 613 cannot be large. B{Bd — ?■ /i^yu^) 
can then distinguish between Bi and B3 oases because cos 613 differs by sign in these two 
oases. We find then destructive interference of Z' contribution with the SM contribution 
in oasis Bi and constructive one in oasis B^ implying the results summarized in Table [6] 
for this \Vub\ scenario. 

We also observe in Table pi that _ can also help by means of its sign to distinguish 

l_l i— I 

between different oases. Fig. 8 (right panel) also shows that in LHS2 the sign of S° 
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Table 6: Oases in the space (sis, 6^) for Mz' = 1 TeV in LHS and two scenarios for 
\Vub\- The enhancement or suppression of B{Bd — )■ fi^fi~) with respect to the SM value 
chooses the oasis uniquely in LHS2 . The sign of S'l^+^- chooses the oasis for both LHSl 
and LHS2. 



LHS1 



LHS2 
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Figure 8: B{Bd fifi) versus S';^+^_ for Mz' = 1 TeV in LHSl (left) and LHS2 (right). 



Bi: yellow, B^: green, B2: magenta, B^: brown. Red point: SM central value. 



is opposite to the sign of the shift in the corresponding branching ratio (except for the 
small oasis i?2(5'2)), which can easily be understood by inspecting the ranges of S13. 
Moreover, the predictions for >S'^+^- are rather precise. This is in particular the case 
for small oases, which in the S^-system cannot be ruled out. In fact in the B2 oasis 
B{Ba — 7- fi^fJ^^) can still be by a factor of three enhanced with respect to its SM value. 
Finally, we observe that S^+^- can be large, although not as large as S^^+^-. 

We next turn to LHSl scenario for \Vub\ which is novel with respect to the analysis 
in 35 . We observe that the phase 613 is lower for big oases than in the case of scenario 



S2 but Si3 is basically the same. We observe that while the sign of S^+^^ can still 
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distinguish between the big oases, 5^ — ?■ /^''"yU" cannot do it as welL This is related to 
the fact that with |\4;,| as low as 0.0031 we have to enhance slightly S^Ks certain 
range of parameters involved. These features are seen in the left panels in Figs. [7] and 

El 

What distinguishes LHSl from LHS2 is the sign of the correlation between 5*1^+^- and 

B{Bd, — >■ A positive 5'^+^- implies enhancement of B{Bd. f^^f^ ) in LHSl 

but suppression in LHS2. Note that this pattern is independent of the sign of Z' ji^ jj,^ 
coupling as this coupling enters both observables. On the other hand the flip of this 
sign would also flip signs in the last two columns in Table [6] and thereby interchange 
colours in Figs. [7] and [8] 



7.1.3 The K Meson System 

In the 331 model, which was governed by S2 scenario, NP effects in rare K decays were 
very small due to suppression of both sdZ' and vvZ' couplings in this model. 

However in a general Z' model the possibility for SI scenario for \Vub\ and enhanced 
values of leptonic couplings with respect to the ones found in the 331 model allow to find 
large NP effects in rare K decays. This allows to find interesting correlations between 
relevant branching ratios that we would like to exhibit here. 



As seen in ( 132 ) the constraints from AF = 2 observables are weaker than in previous 
cases. Yet as seen in Fig. |9]it is possible to identify the allowed oases. These plots have 
the same structure as the plot in Fig. 2 of 34 with the SI and S2 scenario for \Vub\ 



on the left and on the right, respectively. We observe that the small oases are absent 
now as Ek and AM^ are governed respectively by imaginary and real parts of and 
not by their absolute values like in the case of AMg^d- Therefore the solutions with 
very large NP contributions but opposite signs to the SM contributions corresponding 
to small oases in the latter case are not allowed here. 

Due to weaker constraints in the K system the oases are rather large. We have two 
oases in SI: 

Ci{Sl) : 0° < (5i2 < 90°, ^2(51) : 180° < 5^2 < 270° (135) 
and only one oasis in S2: 

Ci{S2): 0°<5i2< 360°. (136) 



As emphasized in 34 of particular interest are the values 



5i2 = n|, 71 = 0,1,2,3 (137) 

for which NP contributions to Ek vanish. As seen in Fig. |9] this is only allowed for 
scenario S2 for which SM agrees well with the data and NP contributions are not 
required. In this scenario 5i2 can even vanish. In scenario SI, in which NP contributions 
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Figure 9: Ranges for AM^ (red region) andex (blue region) (LHSI: left, LHS2: right) 



for Mz' = 1 TeV satisfying the bounds in Eq. (132) 



are required to reproduce the data, Si2 is bounded from below and S12 cannot satisfy 



(137) but for sufficiently large S12 can satisfy it approximately. As at these values of 



S12, the mass difference AMk is non-zero, S12 is bounded from above but due to the 
weak AMx-constraint this is not seen in the plot. 



In 34 an extensive analysis of the interplay between ek and in K ^ txvij in different 



NP scenarios has been performed but the case of tree-level Z' contributions has not be 
discussed. As the latter contributions are much more specific and simpler than the NP 



models discussed in 34 , it will be interesting to see how correlations between Ek and 



K ^ TTw m the eight scenarios in (128) compare with the findings of 34 . 



To this end for the LHS scenarios we find for the quantities defined in [34] 
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(138) 



(139) 



(140) 



Ras=2 2\/r5'SM^z' 

For our choice of A'£'(Z') we find p ~ 0.6 for Mz' = 1 TeV. On the basis of we 
expect for this value of p strict correlation between B{K^ — )■ tx^vv) and B{Kl — )■ tt^uu) 



familiar from the LHT model 32 . It is interesting that p depends only on the size of 



In 34 p was denoted by e. 
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IY^{Z'^ and Mz'- This will have important implications for the study of flavour- 
violating Z couplings considered in Section [91 



In the upper panels of Fig. 10 we show this correlation in LHSl and LHS2 for M^/ 
1 TeV. We observe the following pattern of deviations from the SM expectations: 



There are two branches in both scenarios. The difference between LHSl and LHS2 
originates from required NP contributions in LHSl in order to agree with the data 
on Ek and the fact that in LHSl there are two oases and only one in LHS2. 



The horizontal branch in both plots corresponds to = 0, 2 in (137), for which 



NP contribution to K ^ -nvv is real and vanishes in the case of Ki^ — > vi^vv. 



The second branch corresponds to n = 1,3 in (137), for which NP contribution 



is purely imaginary. It is parallel to the Grossman-Nir (GN) bound 114 that is 
represented by the solid line. 



This pattern agrees with general results of 34 . In fact the structure of plots in the 
Fig. 3 and 4 of the latter paper agrees for p 1 very well with our findings for LHS2 
and LHSl scenarios, respectively. What is striking is the fact that still large deviations 
from the SM predictions are allowed, significantly larger than in the case of rare B 
decays. This is a consequence of the weaker constraint from AS* = 2 processes than 
IS.B = 2 and the fact that rare K decays are stronger suppressed than rare B decays 
within the SM. Yet as we will soon see some of these large values will be ruled out 
through the correlation with Ki — )■ fi~^fi~. 



In the left panel of Fig. 11 we show the correlation between B{K~^ — )■ ir'^uu) and 
B{Kl — >■ fi~^fi~) for LHSl. We note a correlation analogous to the one found in the 
LHT model |32] but due to fewer free parameters in Z' model this correlation depends 
whether oasis Ci or C2 is considered. Very similar correlation is found in LHS2 scenario 
but as here only one very big oases is present only cyan regions appear. We will return 
to the right panel in this figure in the context of RHSl scenario below. 



10 



While at first sight the correlation in Fig. 11 is similar in shape to the one in Fig. 
one should note that Kl — )■ fi'^fj^~ is governed by the real part of the involved master 
function and not imaginary part as was the case of Kl — )■ n^vv. Therefore the horizonal 



line in Fig. 11 corresponds this time to n = 1, 3 in (137), for which NP contribution is 



purely imaginary, while the other branches correspond to = 0,2 in (137), for which 



NP contribution to K ^ ttz/z/ is real and vanishes in the case of tx^vv. 

We observe again that NP effects in both decays can be large and the upper bound on 



B{Kl — )■ /i^/i ) in (68) represented by the horizontal black line can easily be violated. 



The impact of this bound on the results in Fig. [TO] is represented by the black areas 
that violate this bound. 



Combining the information from Figs. 10 and 11 we obtain the following result: 
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Figure 10: B{Kl — )■ ti^vv) versus B{K^ — )► ti^vv) for Mz> = 1 TeV (upper panels, Ci: 
cyan, C2: pink.) and Mz' = 5 TeV (cyan), 10 TeV (blue) and 30 TeV (purple) (lower 
panels) in LHSl (left) and LHS2 (right). Black regions are excluded by the upper bound 
B{Ki — >• fi'^fi~) < 2.5 ■ 10~^. Red point: SM central value. Gray region: experimental 
range of B{K^ n^uu). 



In the case of the dominance of real NP contributions we find in Ci{Sl) for 
Mz' = 1 TeV 

^ TT+z/p) < 16 • 10"". (141) 
In this case Kl ti^uu is SM-like and B{Kl reaches the upper bound 



in (68). On the other hand C2{S1) oasis in this case is excluded through the 



simultaneous consideration of both decays. 

In the case of the dominance of imaginary NP contributions the bound on B{Kl — )■ 
is ineffective and both B{K'^ — tt'^ui?) and B{Kl — tt^uu) can be signifi- 
cantly larger than the SM predictions and B{K~^ — >■ tt^uu) can also be larger than 
its present experimental central value. We also find that for such large values the 
branching ratios are strongly correlated. Inspecting in the LHS2 scenario when 
the branch parallel to the GN bound leaves the grey region corresponding to the 
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Figure 11: B{Kl /x+Zi") versus B{K+ tt+z/P) for Mz' = 1 TeV zn LHSl (left) 
and RHSl (right). Ci: cyan, C2: pink. Red point: SM central value. Gray re- 
gion: experimental range of B{K^ — )■ n^vi') and horizontal black line: upper hound 
ofB{KL^f,+fx-). 



la region in ( 74 ) we find a rough upper bound 

B{Kl vrVz/) < 85 ■ 10-^\ 



(142) 



which is much stronger than the present experimental upper bound in (75). 



We conclude therefore that K — )■ vri/z/ decays provide an important portal to fiavour- 
violating Z' with masses outside the reach of the LHC before its upgrade and even in its 



second phase. In the lower part of Fig. [10] we show therefore how the plots in the upper 
part of this figure would look like for M^/ = 5 TeV, 10 TeV and 30 TeV. We observe 
that even at Mz' = 10 TeV both branching ratios can still differ by much from SM 
predictions and for Mz' < 20 TeV NP effects in these decays, in particular — )■ ir^uu 
should be detectable in the fiavour precision era. For Mz' = 30 TeV and higher scales 
it will be very difficult. 



In the left panel of Fig. 12 we show the correlation between B{Kl — ?■ TT^e~^e ) and 
B{Kl — Tx^n^jj,') that has first been investigated in 89,91,92]. We have shown only 



the results in LIIS2 as our main goal here is to find out whether large enhancements 
of the branching ratios for — ?■ tt^vv and — vr+z/z/ can be affected by these 
decays. To this end we have also assumed constructive interference between SM and 
NP contributions. 



We observe a strong correlation between B{Kl — )■ •n^e'^e ) and B{Kl — )■ 7r°/i ' yU 
similar to the case of LHT 31 and RSc [33j models. Indeed such a correlation is common 



to all models with no scalar operators contributing to the decays in question 89 ,91 ,92 



We also observe that both branching ratios can be in principle enhanced by an order of 
magnitude over the SM values in (86). However, the correlation with — > tt+z/z/ and 
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Figure 12: B{Kl — )• 7r°e+e^) as a function of B{Ki — )• T^^jji^ n^) (left panel) and 
B{Kl — !■ Tr'^e+e^) (upper curve) and B{Kl — 7r°/i+/i^) (lower curve) as functions of 
B{Kl — )■ 7r°z/p) ('ng'/ii panel) in LHS2 for Mz' = ITeV. The red points represent SM 
predictions. 



Kl — -K vv does not allow for such large values. Indeed the experimental upper bound 
on B{K 



rO,„-,^ 



implies in a given scenario an upper bound on B{Ki — )• tt^uu) 
as given in (142) and this in turn implies upper bounds on B{Kl — )■ Tx^e^e^) and 
B{Kl — )■ TT^/i+jU") that are strongly correlated with B{Kl — tt^uu). This correlation 
is evident in Z' scenarios if one compares the expressions for the relevant amplitudes. 



In the right panel of Fig. 12 we show this correlation that has already been found in 
the LHT model [3l] and in the RSc model We note that a large enhancement of 
B{Kl —7- TT^uu) automatically implies significant enhancements of B{Kl — )■ tt^^^^"), 
although the NP effects in Kl — > t^^vv are stronger. This is related to the fact that 
small or moderate NP effects in Kl — )■ are shadowed by the dominant indi- 

rectly CP-violating contribution. However, for large NP contributions when the di- 
rectly CP-violating contribution becomes more important NP effects in B{Kl — )■ 



and B{Kl — )■ 7i^i~^i ) are comparable in size. The correlations in Fig. 
powerful test of the model considered. 

The vertical solid line in the right panel of Fig. 



12 



constitute a 



[12] corresponds to the rough bound in 
(142). It implies upper bounds on B{Kl — > ir^i'^i') that are stronger than the present 
experimental bounds. We indicate these bounds by horizontal and vertical solid lines 
in the left panel of this figure. 

This analysis shows that in LHS scenarios the present upper bounds on B{Kl — )■ 
do not preclude large NP effects found in B{K^ — )■ tt'^ui?) and B{Kl — )■ tt^uu) 
but the bounds on the latter branching ratios have an impact on B{Kl — )■ ir^i^i^). 
This property remains for higher values of Mz'- 
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7.2 The RHSl and RHS2 Scenarios 
7.2.1 First Observations 

We will now investigate Z' scenario with exclusively RH couplings to quarks. We should 
emphasize that such scenario is not artificial as in certain extensions of the SM, the 
corrections to left-handed neutral gauge boson couplings are suppressed due to some 
custodial symmetries. This is for instance the case of the ordinary Z gauge boson in 
Randall- Sundrum scenarios with custodial symmetry in the bulk (RSc). In such a case 



the phenomenology of fiavour violation is dominated by right-handed couplings 33 . 

Now in the RHSl and RHS2 scenarios only RH couplings to quarks are present in Z' 
contributions. As QCD is parity conserving, the hadronic matrix elements for operators 
with RH currents as well as QCD corrections remain unchanged. The expressions 
for AF = 2 observables in RHSl and RHS2 scenarios as well as the corresponding 
constraints have precisely the same structure as in the LHSl and LHS2 cases just 
discussed. Therefore the oases in the space of parameters related to RH currents are 
precisely the same as those given in Tables |5] and [6| except that the parameters Sij 
and bij parametrize now RH and not LH currents. Anticipating this result we have 
not introduced separate description of LH and RH oases. Yet, in the case of AF = 1 
observables several changes are present which allow in principle to distinguish the RHSl 
and RHS2 scenarios from the corresponding LHSl and LHS2, which we just analyzed 
in detail. 

In what follows we will list all changes in the three meson systems one by one. We use 
for the oases the same notation as in the LHS cases. The basic rule for modifications 
of correlations between various observables is as follows: 

• In AF = 2 observables nothing changes as stated above. Therefore we do not 
show any plots for allowed oases. 

• In AF = 1 observables governed by the functions y, that is processes with muons 
in the final state, there is a change of sign of NP contributions in a given oasis. 



See (35) and (36). 



In AF = 1 observables governed by the functions X, that is processes with 
neutrinos in the final state, there is no change of sign of NP contributions in a 
given oasis. The consequences of it are straightforward in the case of — )■ tt^uu 
and Kl — J- tt^uu but in the case of 6 — )■ suu transitions the implications are richer 
as we have four observables to our disposal that are sensitive to the RH currents 



in a different manner. See (78)-(80) and (82). The same comments apply to 
B — >■ K* ix~ and B — )■ Kfi~^fi^ which now receive contributions from primed 
operators Q'g and Q[q. 
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7.2.2 The Meson System 

In Fig. [3] we have shown vs S^^ in the LHSl scenario. This plot is also valid for 

RHSl scenario except that now a negative value of S^+^. chooses scenario Ai, while a 
positive one scenario A^. The size of NP effects is the same in LHSl and RHSl, only 
the oases are interchanged. The same comments apply to large \Vub\ scenario. 

We conclude therefore that on the basis of and S^ff, alone it is not possible to 



distinguish between LHSl and RHSl scenarios because in the RHSl scenario one can 
simply interchange the two big oases or two small oases to obtain the same physical 
results as in LHSl scenario. Therefore let us look at other observables. 

In Fig. [3] (right) we have shown S^^ vs B{Bs — )■ f^'^fJ'') in the LHSl scenario. This plot 
is also valid for the RHSl scenario but again the oases Ai and A3 are interchanged. 
While in the oasis Ai S^^ increases (decreases) uniquely with increasing (decreasing) 
B{Bs — !■ ji^ji'), in the oasis A3, the increase of S^^ implies uniquely a decrease of 

Clearly as in the LHSl scenario this result represents a test of the RHSl scenario but 
if one day we will have precise measurements of S^^ and B{Bs — )• fi^fi') we will 

still not be able to distinguish for instance whether we deal with LHSl scenario in oasis 
Ai or RHSl scenario in oasis ^3. 



Fortunately, as we will see in subsection 7.5, we will be able to make a clear distinction 



between LHS and RHS scenarios by considering model independent bounds from B — )■ 
K* jX^ H~ and B — )■ Kfi'^^~ on the Wilson coefficients of primed operator Q[q. Also, as 



seen in (78)-(80) and (82), b — )■ sz/z/ transitions being sensitive to RH currents will be 
very helpful in this respect. In LHS scenarios the first three observables where affected 
in the same manner and was unaffected. Therefore the comparison of these different 
structures in RHS and LHS scenarios will give us a powerful insight in the Z' couplings. 
As these issues will also be relevant for LR and ALR scenarios, we will discuss b — )■ suu 
observables in the four scenarios in a separate subsection at the end of this section. 



7.2.3 The Ba Meson System 

Similarly to the Bg case the structure of oases is as in Fig. |6] and Table [6j In Fig. [7] 
we have shown S^Ks B{Bd — )• /i'^'/i") for the LHS2 scenario. This plot is also valid 
for RHS2 scenario and B{Bd — )■ yU^/i~) can distinguish between Bi and B3 oases. But 
now the behaviour of B{Bd — t- /U^/i~) in these two oases is interchanged. It is enhanced 
and suppressed with respect to the SM in Bi and B^, respectively. Thus we cannot 
distinguish between LHS2 and RHS2 on the basis of these observables. Clearly the 
study of 6 — )■ and b — > duu transitions could help in this context but they are 

more challenging both theoretically and experimentally. 

Analogous comments apply to RHSl which cannot be distinguished from LHSl on the 
basis of observables considered. 
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7.2.4 The K Meson System 



In Fig. 10 we have shown the correlation between B{K^ — t- n^uu) and B{Kl — > tt^uu) 
in the LHSl and LHS2 scenario. It is evident from ( [73| ) that this plot applies identically 
to RHSl and RHS2 scenarios as well. Thus — )■ n'^uu and — J- ir^uu decays are 
not useful for the search of RH currents as they are sensitive only the vector parts of 
Z' couplings to quarks. As expected, we also find that the results for Kl — > n'^i^fr 
decays are the same as in LHS scenarios. 



However, as known already from different studies, in particular in RSc scenario 33 
the correlation between B{K'^ — )■ vr+i/i/) and B{Kl — )■ fJ''^fJ'~) brings the rescue to this 
problematic as the latter decay is sensitive to the axial-vector couplings. In the right 



panel of Fig. 11 we show this correlation for the RHSl scenario. Indeed the correlations 



in both oases differ from the ones in LHSl. 

We also note that in the case of the dominance of imaginary NP contributions corre- 
sponding to the horizontal line, B{K~^ — > tt'^uu) and B{Kl — )■ ti^vv) can be large. But 
otherwise B{K^ — )■ tt^vv) is suppressed with respect to its SM value and B{Kl tt^vv) 
is SM-like. 

We should emphasize at this point that the impact of Kl — /i^yU^ on — > -k^vv and 
Kl — )■ 7r°z/z/ and its correlation with them depends on the sign of leptonic couplings. 
For a negative the results in LHS and RHS scenarios would be interchanged. 



7.3 The LRSl and LRS2 Scenarios 
7.3.1 First Observations 

If both LH and RH currents are present in NP contributions, the pattern of flavour 
violation can differ from the scenarios considered until now in a profound manner. If 
the LH and RH couplings differ from each other, the number of parameters increases and 
it is harder to get clear cut conclusions without some underlying fundamental theory. 
On the other hand some of the "symmetries" between LHS and RHS scenarios identified 
above are broken and the effect of RH currents in certain cases could in principle be 
better visible. 

Here in order to keep the same number of parameters as in previous scenarios we will 
assume a left-right symmetry in the Z'-couplings to quarks. That is the LH couplings 
are equal in magnitudes and phases to the corresponding RH couplings Ar. In this 
manner we can also keep the same parametrization of couplings as in previous scenarios. 

Before entering the details let us emphasize two new features relative to the cases in 
which either LH or RH couplings in NP contributions were present: 

• NP contributions to AF = 2 observables are dominated now by new LR operators, 
whose contributions are enhanced through renormalization group effects relative 
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to LL and RR operators and in the case of Ek also through chirally enhanced 
hadronic matrix elements. Consequently the oases will differ from the previous 
ones. 

• NP contributions to B^^s — > and Ki^ vanish eliminating in this 
manner 5"^+ _ and B{Bs,d A^^/W") as basic observables in the identification of 
acceptable oases. On the other hand B — )■ K*fi~^fi~ and B — Kfi~^fi~ receive still 
NP contributions and can help in this context. 

• Also NP contributions to decays with neutrinos in the final state, that is — > 
TT+z/z/, Kl — )■ TT^uu and the b — suu transitions are important for testing the 
LRSl and LRS2 scenarios. 

While 5'^+^- cannot help in the identification of the optimal oasis in the LR scenarios 



they are non- vanishing as seen in (59) and (61): 



= -sin(2^sJ. (143) 



While rather small they offer a clean test of LR scenarios. 



7.3.2 The Bs Meson System 

We begin the search for the oases with the Bg system proceeding with input parameters 
as in the previous scenarios. The result of this search for Mz' = 1 TeV is shown in 



Fig. 13, where we show the allowed ranges for (s23, (523)- The red regions correspond to 
the allowed ranges for AM^, while the blue ones to the corresponding ranges for 5*^^. 
The overlap between red and blue regions identifies the oases we were looking for. 

From these plots we extract several oases that are collected in Table [7j The notations are 
as in previous cases but it should be kept in mind that the parameters {§23, S23) describe 
both LH and RH couplings. Note that the entries related to Bg^d Z^"*"/^ are absent 
now as in this scenario there are no NP contributions to these decays. Consequently 
the ATg effects in Bg — )■ fi~^fi are as in the SM: = 1. 

In order to understand the structure of oases in Table [7} that differs from the ones 
found so far, we note that the matrix element of the dominant operator has the 
sign opposite to SM operators. Therefore, this operator naturally suppresses AMg with 
the phase ^23 centered in the ballpark of 0° and 180°, that is shifted down by roughly 
90° relatively to the LHS scenarios. As the matrix element of is larger than that 
of the SM operator in LHS and RHS scenarios, S23 has to be sufficiently small to agree 
with data. 

The crucial role in the Bg meson system in this scenario, in the absence of NP contri- 
butions to Bg d — fi'^fi" decays, is now played by -B — t- K*fi^fi^, B — )■ Kfx^^^ and 
b — 7- suu transitions. We will discuss the latter decays at the end of this section. 
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§23 


^23 




ABf 




0.00059 - 0.00216 


139° - 219° 


± 


T 




0.00628 - 0.00644 


177° - 182° 








0.00059 - 0.00216 


-41° - 39° 


± 


± 




0.00628 - 0.00644 


-3° - 2° 







Table 7: Oases in the space (523,523) for Mz' = ITeV in LRSl. Basically the same 
results are obtained for LRS2 scenario. 



7.3.3 The Ba Meson System 



The structure of oases in this case is given in Fig. [14] and Table [8j As we do not have 
B{Bd — )■ to our disposal and h — )■ dw decays are challenging this system is not 

very useful to provide tests of LRS scenarios without some fundamental theory. 

Due to the sign of the matrix element of the dominant operator in both LRSl and 
LRS2 the mass difference AM^ is naturally suppressed. The requested size of this sup- 
pression together with significant suppression of S^Ks i^i LRS2 and slight enhancement 
of it in LRSl governs the structure of the phases. 
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AMrf a S^K„. LRS1 



AWlrf & S,,K^, LRS2 




0.0000 0.0002 0.00040.0006 0.00030.0010 0.0012 0.0014 



0.0000 0.0002 0.00040.0006 0.00030.0010 0.0012 0.0014 



Figure 14: Ranges for AM^ (red region) and S^Ks (blue region) for Mz' = 1 TeV in 
LRSl (left) and LRS2 (right) satisfying the hounds in Eq. (131). 





Sl3 


^13 


B2{S1) 
BsiSl) 
B,{S1) 


0.00021 - 0.00040 
0.00111 - 0.00117 
0.00021 - 0.00040 
0.00111 - 0.00117 


166° - 195° 
179° - 180° 
-14° - 15° 
-1° - 0° 


Bi{S2) 

B2{S2) 
B3{S2) 
Bi{S2) 


0.00028 - 0.00044 
0.00106 -0.00111 
0.00028 - 0.00044 
0.00106 -0.00111 


38° - 60° 
2° -5° 
218° - 240° 
182° - 185° 



Table 8: Oases in the space (si3,i5i3) for Mz> = 1 TeV and two scenarios for \Vub\ in 
LR scenarios. 



7.3.4 The K Meson System 



In Fig. [15] we show the oases in this system that due to the presence of LR operators 
have a different structure than in previous scenarios. While the shape of the single oasis 
in the LRS2 case is similar to the LHS2, for LRSl the oases are shifted by 90°: 

Ci{Sl) : 90° < 6i2 < 180°, C2{S1) : 270° < 612 < 360°. (144) 
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iSMj( & |sv I . LRB1 



iMjf & lijfl, LRS2 




5. xlO"* 0.00001 0.000015 0.00002 0.000025 0.00003 5. xlO"* 0.00001 0.000015 0.00002 O.OOOOffi 0.00003 

s,2 J, 2 

Figure 15: Ranges for AM k (red region) andex (LRSl: left, LRS2: right) for Mz' = 1 



TeV satisfying the bounds in Eq. (132). 



In Fig. 16 we show the correlation between B{K^ — )■ tt^vv) and B{Kl — )■ tt^z/z/) in 



the LRSl and LRS2 scenarios. Even if this correlation has similar structure to the one 
found in previous scenarios, there are differences: 



Concentrating first on the Mz' = 1 TeV case we observe that the branches in both 
scenarios are much thiner than in the LHSl and LHS2 cases which originates in 
the fact that NP contributions to Ek have to be kept under control in the presence 
of LR operators. 



In the LRSl case the new structure of phases in (144) implies definite predictions 



for both branching ratios in Ci(5'l) and (72(5*1) respectively (M^/ = ITeV): 

3-10"" < B{K+ n^uu) < 12-10"^\ 2-10"^^ < B{Kl -> vrVi^) < 14-10"", 

(145) 

8 ■ 10"" < B{K-^ TT+z/z/) < 16 ■ 10"", < B{Kl ir^yv) < 3 ■ 10"". 

(146) 

In LRS2 the effects are slightly larger than in LRSl and in fact on the horizontal 



line B{K^ — t- -k^vv) can be larger than in LHS2 in Fig. 10 as in LRS2 the 
Kl — /i^/i" bound is ineffective. 

However, otherwise the effects in LRSl and LRS2 scenarios are much smaller than 
in LHSl and LHS2 scenarios in accordance with the correlations between AF = 1 



and AF = 2 transitions derived in Subsection 5.2 This is clearly seen in lower 
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LRS1 



LRS2 




BtK^'^ff' V V) [ID""] 



5 10 15 20 25 30 
B(K"^ff" 1' V] [ID"''] 



LRS1. Cyan: 5TeV, Blue; lOTeV. Purple; 30TeV LRS2. Cyan; 5TeV. Blue; IGTeV. Purple; 30TeV 

7 l J 7r 




BcK""^ff" 1^ 1') [10""] 



7 8 9 

B(K"^ff" imO 110""] 



10 



Figure 16: B{Kl — J- n^uu) versus B{K^ — t- ti^uu) for Mz' = 1 TeV in LRSl and 
LRS2 (upper panels, Ci: cyan, C2: pink.) and Mz' = 5 TeV (cyan), 10 TeV (blue) 
and 30 TeV (purple) in LRSl and LRS2. Red point: SM central value. Gray region: 
experimental range of B{K^ — n^uu). 



panels in Fig. 16, where we show the results for higher values of Mz'- While for 
Mz' = 5 TeV NP effects in both branching ratio can still be distinguished in the 
future from SM values, for higher masses of Mz' this will be very difficult. Note 
that the SM values of B{Kl — )■ tt^uu) in LRSl are visibly smaller than in LRS2 
which is a clear consequence of the change of \Vub\- 

This discussion shows that if Mz' > 5 TeV and both branching ratios will be found 
significantly larger than SM values the LHSl and LHS2 scenarios will be favoured over 
LRSl and LRS2. 

These results have been obtained under the assumptions of the exact LR symmetry. 
However, as already shown in |115 and also discussed in 34 if the phase structure of 



LL and RR and LR contributions in AS* = 2 transitions is not related to each other 
the correlation between B{K^ — )■ ti^vv) and B{Kl — )■ tc^uu) can change profoundly. 
This is for instance seen in Fig. 3 of |115| , where for such a general scenario one can 
find a decrease of B{Ki — )■ -k^uj?) with increasing B{K~^ — > tt^uu), the property which 
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LR32 



LRS2 




B{Ki_^iz' ji' ii' ) [10""l 



Figure 17: B{Kl — vr^e^e^) as a function of B{Kl — 7r°yU+/i^) (left panel) and 
B{Kl — > 7r°e+e^) (upper curve) and B{Kl — )■ 7r°/i"'"/i~) (lower curve) as functions of 
B{Kl — )• n^uu) (right panel) in LRS2 at Mz' = ITeV. The red points represent SM 
predictions. 



is absent in the results presented sofar. Indeed in this case the leading NP contribution 
in AS* = 2 amplitudes is not proportional to the square of NP contribution to AS* = 1 
amplitudes implying a different correlation. As such a study can be more efficiently 
performed in a concrete model, we leave it for the future. 

Finally we discuss Kl 7r°£"'"£~ decays in this scenario. In Fig. 17 we show the results 



corresponding to those in Figs. 12 obtained in LHS2. We observe that NP effects in 
Kl — n^i'^i' in this scenario can be large but much smaller than in LHS2 in accordance 
with the correlations derived in Subsection 5.2, We have also checked that the present 



bounds on these decays do not remove sizable NP effects found in LRS2 scenario for 
— !■ TT'^uu and Kl — > n^uu decays. Similar results are found for LRSl. 



7.4 The ALRSl and ALRS2 Scenarios 



We include this case as well because it has not been discussed in the literature but it is 
an interesting NP scenario for the following reasons: 

• NP contributions to AF = 2 observables are dominated as in LRS scenarios by 
new LR operators but as the sign of LR interference is flipped some differences 
arise. 

• NP contributions to Bd^s enter again with full power. Therefore these 
decays together with S^^^^^ offer as in the LHS and RHS scenarios some help in 
the identification of acceptable oases. 

• The phase structure of the oases is as in LHS scenario but due to enhanced 
hadronic matrix elements of LR operators the mixing parameters Sij are decreased 
typically by a factor of 3.5. 
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• NP contributions to — > -n^vv^ Kl — )■ n^vv and Kl — )■ ^l^t^^ vanish in this 
scenario but the h — ?■ svv transitions can still offer important tests. 

In view of this simple structure of modifications with respect to LHS scenario, all plots 
with the exception of decays with neutrinos in the final state have the same structure 



as LH scenarios but NP effects are smaller. This is also seen by inspecting (107), (108) 



(114) and (115). Therefore we will not show these plots. However when we will move to 
consider flavour-violating Z couplings this suppression of NP effects in ALRS relative 
to LHS will turn out to be fortunate and ALRS will be doing better than LHS. 

Concerning NP contributions to K ^ vrz/z/ decays this scenario could turn out one day 
to be interesting if the data on observables in Bg and Bd systems will show the presence 
of RH currents but negligible NP effects m K ^ -nvv. On the other hand as we will 
discuss soon h — svv transitions can help to distinguish this scenario from the previous 



ones. 



7.5 Implications of 6 — > st^i Constraints 



Presently the NP effects found by us are consistent with the experimental data on 
f^~^f^~y although a range of values above the SM estimate of B{Bs — )■ fi^fi") 
already slightly violates the existing upper bound. However, also the data on i? — t- 
Xsi~^i~ , B — )■ K*l^l~ and B — )■ K£^£^ improved recently by much and it is of interest 
to see whether they have an impact on our results. A very extensive model independent 
analysis of the impact of the recent LHCb data on the Wilson coefficients Cg''' and C^''^ 



10 



has been performed in 36 and we can use these results in our case. 



3.2.3 



the Wilson coefficients Cg''' depend on Ay^{Z') couplings 
They can be chosen to satify the constraints in 



As seen in Subsection 
which did not enter our analysis, 
question. Therefore we will only check whether for the ranges of parameters considered 
by us the resulting coefficients CfJ satisfy the model independent bounds in 



36 



As 

these coefficients are scale independent we can use the formulae in Subsection |3.2.3| and 
compare the resulting coefficients with those in the latter paper. The allowed 2a ranges 
of C|q are shown in Figs. 1 and 2 of j36|. They are given approximately as follows: 

(147a) 



-2<^{C[q)<0, -2.5 < 3(C;o) < 2.5, 



- 0.8 < 3ft(Cfo^) < 1.8 , -3 < 53(Cio) < 3 



(147b) 



Especially, the new data on B ^ K*fi^fi allow only for negative values of the real 
part of C[q 

^iC[,) < (148) 
and this has an impact on our results in RH and LR scenarios. 
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even 



We find that the NP effects in C\q are within the la bounds presented in 36 
if with improved data one will be able to remove certain range of parameters. Some 
footprints of this can already been seen in Figs. [I and H} 



In Fig. 18 we show the results on C{q in RH and LR scenarios without the constraint 
in Eq. (148). We observe that imposing it removes roughly half of each of oases A\ and 
A3 in the RH case and the oasis Ax (blue) in the case of LR. We will now investigate 
the implications of this finding on the comparison between LHS and RHS scenarios. 

The most interesting impact of this constraint at present is on the correlation of S^^ and 
B{Bs Ai^A* ) that for the LHS is given in Fig. |3} As we discussed previously the same 
result is obtained in the RHS with the two oases A\ (blue) and ^3 (purple) interchanged. 
However taking the constraint ( 148 ) into account results in a modified correlation within 

The black areas are excluded and in the RH 
relative to the SM is excluded. This 



19 



RH scenario that we show in Fig. 
scenario an enhancement of BiBg 
feature distinguishes LH and RH scenarios. For E{Bs A*^/^ ) below its SM value, the 
measurement of these two observables cannot distinguish between LH and RH scenarios 
as one can always move to the other oasis to obtain the same result. 



7.6 h — )- sw Observables in Different Scenarios. 

In view of the important role of these transitions in the search for RH currents we 



devote to them a separate subsection. We begin with the e — 77 plane proposed in 83 



In Fig. 20 we show the results for all four scenarios considered by us. Indeed they 
can be clearly distinguished in this plane. Indeed a future determination of e and t] 
will tell us whether the nature chooses one of the scenario considered by us or a linear 
combination of them. 

With four observables and four scenarios for Z'-couplings there is a multitude of results 
for specific observables one could present at this stage. Here we present only some of 
them that we consider most interesting. 
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RHS1 




2 3 4 5 



Figure 19: S^^ and B{Bs — )■ ) for Mz' = 1 TeV in RHSl. Ai: blue, A^: purple, 
A2: red, A/^: gray. Black: excluded due to '^{C'^^ > 0. Gray region: exp la range 
B{Bs At+/u") = {2.9t\i) ■ 10"^- Red point: SM central value. 
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Figure 20: rj versus e for scenario LHSl, RHSl, LRSl and ALRSl. 



On the left in Fig. 21 we show B{B — )■ Kvv) versus B{Bs — ?■ /i^/i^) in LHSl and 
RHSl scenarios. Even without the ( 148[ ) constraint which ehminates the black region 
in the RH scenario, there is a clear distinction between these two scenarios so that the 
measurement of these two observables can uniquely tell us whether we deal with LHS 
or RHS case. Imposing (148) we find that in the RH scenario B{B — )■ Kvv) can only 
be enhanced and B{Bs -> /i^/i^) suppressed. In the case of B{B — )• K^uu) as seen on 
the right in the same figure its branching ratio can only be suppressed relative to the 
SM in RHS but otherwise the distinction between LHS and RHS is not as pronounced 
as for B{B ^ Kvv). 



In Fig. 22 we show the correlations between B{B — > Kuv) and S"^,^ (top) and B{B — 
K*i>u) and S^,^ (down) in LHS and RH scenarios with the black regions excluded by 
the constraint in (148). We note in particular that in the RHS the measurement of S^p^, 
if different from the SM value, will uniquely determine the allowed oasis. 
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LHS1 (blue) & RHS1 (brown) LHS1 (lue) & RHS1 (brown) 




Figure 21: E{B — )■ Kvv) versus B(Bs — )■ /i"*"/^") (left) and B{B — t- K*vi)) versus 
B{Bs — )■ (right) for Mz' = 1 TeV in LHSl (blue for both oases Ai^^) and RHSl 

(brown for both oases Ai^^)), black points excluded by b ^ si^i~ . Gray region: exp la 
range B{B, fi+fi-) = {2.9t\i) ■ 10"^- -Red point: SM central value. 



LHSl 



RHS1 




Figure 22: B{B — )• Kuu) versus S^^ (top) and B{B K*uu) versus S^^ (down) for 
Mz' = 1 TeV in LHSl (left) and RHSl (right). A^: blue, Aj,: purple, A2: red, A^: 
gray. Black points excluded by b ^ si~^i~ . Red point: SM central value. 
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LHS1 (blue) & RHS1 (brown) LHS1 ( bfue) & LRS1 (purple) 




B( B^KPv)/B(B^KriOEM B( B ^Ku i')/B( B ^Kr r) 

Figure 23: B{B K*uu) versus B{B Kuu) for Mz' = 1 TeV in LHSl (blue for 
both oases ^1,3^, RHSl (brown for both oases ^1,3^ and LRSl (purple for both oases 
black points excluded by b ^ si^i~ . Red point: SM central value. 



Finally, in Fig. 23 we show B{B — )■ K*vv) versus B{B — )• Kuu) in RH and LR scenarios 
compared to LH scenario with black regions excluded by the constraint in (148). Also 



these plots will test the presence of RH currents as emphasized in 115 



The observable Fi being dependent only on rj can in principle serve to identify the 
presence of RH currents. However in our case, as seen in Fig. 20, I77I < 0.16, and NP 



effects in Fl turn out to be small as can be deduced from Fig. 4 in 83 



8 The U(2f Limit 



We have investigated how the parameter space is further constrained when the flavour 
f/(2)3 symmetry 23 29 is imposed on the Z' couplings. As pointed out in 30 in this 
case (fiBa = fBs which in turn implies not only the correlation between CP asymmetries 
S^Ks and S^cf) but also a triple S^Ks ~ S^c/, — \ Vub\ correlation. 

As in [30j we will only consider the case of [7(2)^ broken by the minimal set of spurious: 

ij 



the MU{2)^ case. Then only the LHSl and LHS2 are involved. We find then that s 



and Sij are constrained at the fundamental level as follows. 



In the K system we have 



S12 = a|\4rf||Ft 



ts I ) 



6 



12 



(149) 



where a > and real. Thus NP effects in ek, — ?■ tt+z/i/, Kl — )■ ti^uu and Kl — > fj^~^fi~ 
are described for fixed leptonic couplings by a single real and positive definite parameter. 
Once this parameter is fixed through one of these observables, the others are uniquely 
predicted. Note that this is even more predictive than CMFV in which except for 
common CKM couplings there is no relation between AF = 2 and AF = 1 transitions 
unless some ratios are constructed 116 
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The observables in Bd and Bg systems are correlated with each other due to the relations: 

5l3 S23 



\Vtd\ \Vts\ 



Sl3 - ^23 = /3 - Ps. (150) 



Thus, once the allowed oases in the Bd system are fixed, the oases in Bg system are 
determined. Moreover, all observables in both systems are described by only one real 
positive parameter and one phase, e.g. (5235^23)- 

In Fig. 24 we combine Figs. ^ and fusing the U{2Y symmetry relations in (150). In the 



f/(2)3 limit the small oases are eliminated and the big oases get smaller. This decrease 
turns out to be not very pronounced in the case of (si3, ^13) oases as they were already 
small as seen in Fig. |6] but has a significant impact on {§23, 623) oases which where 
much larger as seen in Fig. [2j Moreover the fact that the results in the Bd system 
depend on whether LHSl or LHS2 is considered is now transfered through the relations 
in ( 150[ ) into the B^ system. This is clearly seen in Fig. 24, in particular the final oases 



(magenta) in LHS2 are smaller than in LHSl due to the required shift of S^Ks- 
This change of allowed oases in the Bg system has a profound impact on the correlation 



between S^^f, and B{Bs — > /i"^;U~). We show this in Fig. 25 that should be compared with 
the corresponding correlation in Fig. |3j We observe that already the sign of 5*^0 will 
decide whether LHSl or LHS2 is favoured. Moreover if B{Bs — )■ fj.~^fi~) will turned out 
to be suppressed relatively to the SM then only one oasis will survive in each scenario. 
Comparison with future precise values of \Vub\ will confirm or rule out this scenario of 
NP. These correlations are particular examples of the correlations in MU{2Y models 



pointed out in 30 . What is new here is that in a specific model considered by us the 



\Vub\ — S^^ correlation has now also implications for B{Bs — )■ /i^/x 



We also note that in this case 30 



S;^,- = = M'^^Y - 2y.new), (151) 

where 

^y = 4 = ^y, ^new = = VB,- (152) 

Moreover, as the CMFV relations for B{Bs,d — > /^^Ai ), also apply in this case the result 



in (65) allows to find (30j 

B{Bd /i+/i") = (1.0l°J) X 10"^°, (CMFV, MU(2)='). (153) 



Finally we remark that within the 331 model, analyzed by us in 35 , the imposition of 
U (2)3 symmetry implies the relations: 

Si3 = b\Vtd\, §23 = b\Vts\, 6,-62 = P-l3s + 7r, (154) 

where b > and real. Consequently NP in all three systems is described by only one 
real positive definite parameter and one phase. 
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AM, i & AMrf a S^Ks> LHS1 U(2)-' 



AM, & S^^ a AMrf S S^K^, LHS2 U(2)~ 




0.000 



Figure 24: Ranges for AM^ (red region), S^^ (blue region), AM^ (green region) and 
Si/jKs (yellow region) for Mz' = 1 TeV in LHSl (left) and LHS2 (right) in the [7(2)^ 



limit satisfying the hounds in Eq. (130) and (131). The overlap region of all four 
regions is shown in magenta. 
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Figure 25: S^^ versus B{Bs fi+fi') for Mz' = 1 TeV m LHSl (left) and LHS2 
(right) in the U{2)^ limit. Blue region corresponds to the lower magenta oases in Fig. 24 
(former Ai) and the purple region corresponds to the upper magenta oases (former A3). 
Gray region: exp la range B{Bs — )■ = (2.9]'^j;'^) ■ 10~^. 
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9 Flavour Violating SM Z Boson 



9.1 Preliminaries 



We will now turn our attention to flavour violating Z-couplings that can be generated 
in the presence of other neutral gauge bosons and or new heavy vectorial fermions 



with +2/3 and —1/3 electric charges. RSc is an explicit model of this type 33,117 



See also 18 . Recently, an extensive analysis of flavour violation in the presence of a 



vectorial +2/3 quark has been presented in 118 , where references to previous literature 



can be found. In the case considered by us, new quarks with —1/3 charges are essential 
for generating flavour violating couplings to SM down-quarks but the presence of heavy 
quarks with +2/3 charges could be relevant for charm physics. Moreover, such heavy 
fermions could contribute to rare K and B decays through loop diagrams. In what 
follows we will assume that these loop contributions can be neglected in comparison 
with the tree-level effects discussed by us. Of course in a concrete model one has to 
check whether this assumption is justifled. 

The strategy and formalism developed in the previous sections can be used in a straight- 
forward manner for the case of Z flavour-violating couplings to quarks. In this case we 
have 

Mz = 91.2 GeV, Af (Z) = Af{Z) = 0.372, A>f{Z) = -0.028 (155) 
The implications of these changes are as follows: 



• The decrease of the neutral gauge boson mass by an order of magnitude relatively 
to the nominal value Mz' = 1 TeV used by us decreases the couplings Sij by the 
same amount without any impact on the phases 6ij when the constraints from 
AF = 2 processes are imposed. 

• As already noticed in ^5] and discussed at the beginning of our paper once the 
parameters Sij are constrained through AF = 2 observables the decrease of neutral 
gauge boson mass enhances NP effects in rare K and B decays. This follows from 
the structure of tree-level contributions to FCNC processes and is not generally 
the case when NP contributions are governed by penguin and box diagrams. The 
formulae in Section [5] exhibit this feature transparently. 

• The latter fact implies that already the present experimental bounds on B{K^ — )■ 
TT+z/z/) and B{Bs,d /^^yU") as well as the data on i? -> X^f+f", B K*i+i~ 
and B — t- Ki~^i~ decays become more powerful than the AF = 2 transitions in 
constraining flavour violating couplings of Z so that effects in AF = 2 processes 
cannot be as large as in Z' case. 



We will now investigate how the flavour-violating Z couplings perform in the three 
meson systems. 
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9.2 The Bs Meson System 

We used first only the AM^ and S^^ constraints finding that small oases are excluded 
by the data on Bg Z^"*"/^"- However, also in big oases B{Bs — ^ A*"''/^") is always 
larger than its SM value and mostly above the data except in LRS case where NP 
contributions vanish. In Fig. [26] (in the case of the LHSl scenario) we show the region 
allowed by AM^ and S^^ contraint together with the yellow region allowed by B{Bs — )■ 
/i^/i") = (2.9;!;}'^) ■ 10~^. We observe no overlap between these regions. A small orange 
region is still left when 1o range for BiB^ A^"*"/^") is considered. Thus at first sight 
one could conclude that in contrast to the Z' scenarios B{Bs A^^/^ ) could reach 
values 2cr away from its experimental value. 

However, when the values of C^^ and CJg considered the situation gets worse. The 
2cr bounds on these coefficients from [36] (see Eq. ( 147[ )) are violated in LHS and RHS 



cases and for the LRS scenario only small room is left. Indeed we find 

|^>(Cio)| > 3.5 (LHSl), |$>(Cio)'| > 3.5 (RHSl) (156) 

and 

|5R(CfoP)| > 1.5 (LRSl), |3?(Cio)'| > 1.5 (LRSl). (157) 

Therefore, the main message from this exercise is that when the above constraints are 
taken into account it is very difficult to suppress AM^ sufficiently in LHS, LRS and RHS 
scenarios without violating the constraints from h — ?• sfi~^^~ transitions. We conclude 
therefore that this NP scenario appears to be strongly disfavoured even if not fully ruled 
out because of assumed small hadronic uncertainties. 

In ALR scenario we find 

|3ft(Cfo^)| > 1-5 (ALRSl), |3(Cio)'| > 1.0 (ALRSl). (158) 

While this scenario is in a slightly better shape most of the allowed space is ruled out 
as well. 

We have also calculated Cg and Cg coefficients. Due to the smallness of A^* in the SM, 
the present 2a bounds are satisfied in all scenarios. Yet, in view of the results for Cio 
it does not look that we should expect much from flavour-violating Z couplings in Bs 
system and consequently we will not consider b — )■ suu transitions. Similar conclusions 



have been reached in 36 37 . 



9.3 The Bd Meson System 



In the Bd system using the same constraints as before we find the allowed oases shown 
in Fig. 27 While the magenta regions in this plot are also allowed by the upper bound 
on B{Bfi — i- |U+/i~), Fig. 28 shows that the latter bound has already and impact on 
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AM; & S..^ & Bj^/iV. LHS1 
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0.0010 



0.0015 



Figure 26: Ranges for AMg (red region), S^^ (blue region) and Bg — )■ /^'''/x^ (yellow) in 
LHSl satisfying the hounds in Eq. (13()) and B{Bs — )■ jJ^'^/J'') in la range [1.8,4.3] ■ 10~^ 
(yellow) and [0.7,5.7] • 10~^ (orange). 



the LHSl and LHS2 scenarios. This should be compared with the Z' case in Fig. [7| 
where NP effects where much smaller. This time we show also the results in ALRSl 
and ALRS2 scenarios in which NP effects are smaller than in LHSl and LHS2 scenarios. 
We only show the results for big oases as small oases imply typicaly values for B{Bd 
//"•"/i") of 8 • 10^^ and 4 ■ 10~^ for LHS and ALRS, respectively and are ruled out. 

We observe that B{Bd — )■ /J^'^IJ'') can strongly be enhanced in all shown scenarios so 
that with improved bound on this branching ratio LHSl and LHS2 scenarios could be 
put into difficulties, while in ALRSl and ALRS2 one could easier satisfy these bounds. 
If such a situation really took place and NP effects would be observed in this decay, this 
would mean that both LH and RH Z'-couplings in the Bd system would be required 
but with opposite sign. 



Fig. 29 shows that not only B{Bd ) but also the CP-asymmetry S'^+^^ can 



deviate significantly from SM expectation. 



9.4 The K Meson System 

The effects of flavour violating Z couplings in tx^vv and Kl — > ix^vv can 

be very large in LHS, RHS and LRS but they can be bounded by the upper bound 
on Kl — )■ except for the LR scenarios and the case of purely imaginary NP 

contributions in all these scenarios where this bound is ineffective. 
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Figure 27: Ranges for AMd (red region), S^Ks (blue region) and Bd — > ^ (yellow) in 
LHSl (left) and LHS2 (right) satisfying the hounds in Eq. (131) and B{Bd — ?■ /i^/i^) < 
9.4-10-1°. 



We begin therefore with LRSl and LRS2 scenarios and show in Fig. 30 B{Kl 
vs B{K^ — TT^vv). Indeed the NP effects can be much larger than in the Z' case shown 
in Fig. [161 



In Fig. 3]^ we show analogous result for LHSl and LHS2 case imposing the Kl — )■ yU^/i^ 
constraint. We observe that in both scenarios only the branch unaffected by the Kl — > 
constraint survives. This is C2{S1) and Ci{S2). In Ci{Sl) B{K^ — t- -n^uu) > 
3.5 ■ 10-1° and this case is ruled out. In Fig. p2^we show the corresponding results for 
RHSl and RHS2 where the — )■ ^~ constraint has a different impact than in LHS 
cases. Also here NP effects can be very large. 

rO/ 



Finally we discuss Ki — )■ vr^f+f" decays. In Figs 
corresponding to the ones found for Z' LHS and LRS scenarios 
SI scenarios turn out to be very similar we only show the results for S2 scenarios. 



33 and 34 we show the results 
As the results for 



The meaning of the curves is the same as in the case of Z' results in Figs. 12 and [T7| 
We observe that NP effects in Kl — )■ decays can be very large but they are 

bounded by the upper bound on B{Kl — )■ ti^vu) which follows from the present bound 
on B{K~^ — )■ n^uu). We note that in LHS scenarios the upper bound on B{Kl tt^uu) 
practically coincides with the GN bound and amounts to 



B{Kl -> vrVz/) < 115 ■ 10" 



-11 



(159) 



It is slightly weaker in LRS scenarios. In any case the the present upper bounds on 



BiK^ 



— > tt'- 



do not preclude large NP effects found in B{K^ — > n^uu) and 
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Figure 28: S^Ks versus B{Bd — )• ) m LHSl, LHS2 (upper row) and ALRSl, 

ALRS2 (lower row) . Bi: yellow, B^: green. Red point: SM central value. Gray region: 
excluded by B{Bd fi+fi') < 9.4 ■ 10"^°. 



B{Kl ttVz/). 



9.5 Comments on e'/e 

The large NP effects in — )■ tt'^uj? and Kl — t- tt^uj? both through Z' and Z tree level 
exchanges belong clearly to highlights of our paper. Yet we would like to emphasize 
at this point that in principle these large effects could be eliminated by the ratio e'/e 
if the relevant hadronic matrix elements were precisely known. As already pointed 
out in 



44 



45 there is a strong correlation between — )■ n^vu and Ki — )■ 



and e'/e because electroweak penguin contributions that are relevant for e'/e govern 
K'^ — i- vr+z/z/ and Kl —¥ tt^ui/ even if operators are different. The strongest correlation 
is between B{Kl — )■ tt^uu) and e'/e because they are both CP- violating. However, if 
B{Kl — )■ 7r°z/i/) is bounded by e'/e then automatically B{K~^ — )• n'^uu) is bounded on 
the branch parallel to the GN bound on which their ratio is approximately constant. On 
the second branch B{K~^ — t- tt^uu) is less affected but there the bound from Kl — ?■ yU"''/^" 
plays a role unless we work in LRS scenario. 
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LHS1 LHS2 




Figure 29: B{Bd — )• ) versus S^+^- LHSl, LHS2 (upper row) and ALRSl, 

ALRS2 (lower row). Bi: yellow, B^: green. Red point: SM central value. Gray region: 
excluded by B{Bd < 9.4 • 10-^°. 
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Figure 31: E{Ki^ vrVz/) versus B{K+ Ti+vv)m LHSl (left) and LHS2 (right). Ci: 
cyan, C2: pink. Black points excluded due to B{Ki — )■ ii^ii") < 2.5 ■ 10~^ constraint. 
Red point: SM central value. Gray region: experimental range of B{K^ ti^uu). 
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Figure 32: B{Kl -k^vv) versus B{K^ 
cyan, C2: pink. Black points excluded due 
Red point: SM central value. Gray region: 



RHS2 




20 40 60 80 
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K+uu)in RHSl (left) and RHS2 (right). Ci: 
to B{Kl — i- fi^fi^) < 2.5 • 10^^ constraint, 
experimental range of B{K^ — )■ -k^vv). 



Unfortunately, in spite of the recent progress on the calculation of the hadronic matrix 
elements relevant for e' /e 119 121 , the hadronic uncertainties in e' /e are still too large 
for reaching a clear cut conclusion on the impact of this ratio on rare K decays. An 

If 



46 



analysis of e' /e in the LHT model demonstrates this problem in explicit terms 
one uses hadronic matrix elements of QCD and electroweak penguin operators obtained 
in the large N approach, {e' /e)Q,yi is in the ballpark of the experimental data and sizable 
departures of B{Kl — )■ tt^uu) from its SM value are not allowed. — t- tt^uu being 
CP conserving and consequently not as strongly correlated with e' /e as B{Kl — )■ n^vv) 
could still be enhanced by 50% in LHT. On the other hand if hadronic matrix elements in 
question differ significantly from their large N values, (£'/£)sm disagrees with experiment 
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Figure 33: B{Kl — )• vr°e+e^) as a function of B{Kl — )• 7r°/i+/i^) ('/e/i panel) and 
B{Kl — !■ TT^e+e^) (upper curve) and B{Kl — )■ 7r°/i+/i~) (lower curve) as functions of 
BBt^Kl — n^uu) (right panel) in LHS2. The red points represent SM predictions. 
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Figure 34: B{Ki — )• 7r°e+e^) as a function of B{Ki — > 7r°yU+/i~) ('/e/^ panel) and 
B{Kl — )• 7r°e+e~) (upper curve) and B{Kl — )■ n^jX^jj,^) (lower curve) as functions of 
BBr{KL — )■ n^uu) (right panel) in LRS2. The red points represent SM predictions. 



and much more room for enhancements of rare K decay branching ratios through NP 
contributions is available. 



10 Summary and Conclusions 

In this paper we exhibited the pattern of flavour violation in models in which NP effects 
are dominated by tree-level Z' exchanges under the assumption that the theoretical and 
experimental errors on various input parameters will decrease with time. In particular 
we have identifled a number of correlations between AF = 2 and AF = 1 processes 
that will enable in due time to test this NP scenario. Our detailed analysis of these 
correlations in Section [7] shows that in the three meson systems considered signiflcant 
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Bd,s systems: observables 


scenarios 


Figure 


comments 


AMg & S^(j, oases 


LHS/RHS Z' 
LRS Z' 

LHS C/(2)3 Z' 
LHS Z 




2 
13 
24 
26 




oases {Ai, A2,A:i, An) found 
oases (Ai, A2,A:j, An) found 
two common oases for Bs and 
no oases found 


5^0 vs. Bs /U+;U" 


LHS Z' 

RHS Z' 
LHS [7(2)3 Z' 




H 

19 
25 




excludes {A2, A4) 
(anti)correlation found in A3 (Ai) 
(anti)correlation found in Ai (A3) 




LHS Z' 




3 






B —7- Xgvv vs. -Bs — >• IJ'~^fJ'^ 


LHS Z' 




4 




correlation found in A3, Ai 


B — )• Xsi/j/ vs. 


LHS Z' 




5 




(anti)correlation found in A3 (Ai) 


K^^^uu vs. -Bs — > 


LHS/RHS Z' 




21 




B -> KPi^ : {anti)correlation in LHS(RHS) 
B — >■ K*vv : correlation in LHS and RHS 


-B ^ i^(*)i>z^ vs. 


LHS/RHS Z' 




22 




LHS: (anti)corrclation in A3 (Ai) 
RHS: K case: above SM, 
(anti)correlation in A3 (Ai); 
opposite in K* case: below SM 


-B — ^ i^*Pz^ vs. -B — ^ Kj^v 


LHS/RHS/LRS Z' 




23 




(anti)correlation in LHS and LRS (RHS) 


/m(C(o) vs. fe(C(o) 


RHS/LRS Z' 




18 




RHS: /m(CJo) > (< 0) in A3(Ai) 
Re(c;o) > 0(< 0) in Ai(A3) 


?7 vs. e (6 — )■ sPz^) 


LHS/RHS/LRS/ALRS Z' 




20 




»j = in LHS; no dependence on e in RHS 
(anti) correlation in ALRS(LRS) 


AMd & S^Ks oases 


LHS/RHS Z' 
LRS Z' 

LHS ?7(2)3 Z' 
LHS Z 




6 

14 
24 
27 




oases (Bi, B2,B3, B4) found 
oases (Bi, B2,B3, B4) found 
two common oases for S3 and Bj 
only two oases found 


•S'^i^s vs. -Bd 


LHS Z' 
LHS/ALRS Z 




H 

28 




B{Ba -> M+At~) above SM in B2, B4, 
and in Bi (LHSl) or in B3 (LHS2) 
Bii — > fJ,'^ /J,^ always above SM; 


5, ^ vs. 
i?, ^ vs. 


LHS Z' 
LHS/ALRS Z 




8 

29 




^U,,- >0(<0) in Bi,B4{B2,B3) 



Table 9: Overview of correlation plots in and Bg sectors. 



foot prints of Z' will be seen provided M^/ < 5 TeV. But only in the K system these 
effects can be detected if Mz' is larger and outside the LHC reach. 

In view of these findings we have investigated whether larger effects could be found 
if the SM Z boson couplings were flavour violating. Indeed, in this case the stronger 
constraints come from AF = 1 and not AF = 2 processes. We find then that imposition 
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K system: observables 


scenarios 


Fi 


gure 


comments 


AMx & £k oases 


LHS/RHS Z' 




9 




oases (Ci,C2) found in 81; only Ci in S2 




LRS Z' 




15 




oases (Ci,C2) found in SI; only Ci in 82 


Kl — )• -K^vu vs. K'^ — )• 'K'^VV 


LHS Z' 




10 




two branch structure 




LRS Z' 




16 




two branch structure 




LHS Z 




31 




only one branch allowed 




LRS Z 




30 




two branch structure 




RHS Z 




32 




two branch structure 


Kl — )• fi'^fi~ vs. -fC+ — vr+z^z^ 


LHS/RHS Z' 




11 




(anti)correlation in LHS (RHS) 


-fCi — )• vr'^e^e" vs. Kl — ?• ■k^ijl'^ j-r 


LHS Z' 




12 




correlation 




LRS Z' 




17 




correlation 




LHS Z 




33 




correlation 




LRS Z 




34 




correlation 


Kl 7r°^+^- vs. ^ 7r°z?i/ 


LHS Z' 




12 




correlation 




LRS Z' 




17 




correlation 




LHS Z 




33 




correlation 




LRS Z 




34 




correlation 



Table 10: Overview of correlation plots in the K sector. 



of the present constraints from h — > s/i^/i^ transitions makes NP effects in AF = 2 
transitions in the Bg system very small precluding also large NP effects in 6 — )■ svv 
transitions. 

The situation is different in the Bd system and in particular in the K system where 
effects from flavour- violating Z couplings are still allowed to be large. 

Our results are summarized in a number of plots that have been obtained in various 
scenarios for the Z' couplings and for inclusive and exclusive values of \Vub\- Also a 
number of plots have been shown for the case of flavour- violating Z- couplings. The 
overview of all correlations found by us and of the related figures is given in Tables [9] 



and [TOj There we collect a short description of all those plots, in particular stressing 
whether correlations or anticorrelations are found among the various observables. We 
list here only few highlights: 



• For each scenario we have identified allowed oases in the parameter space of 
the model. In each oasis particular structure of correlations between various 
observables will in the future either favour or exclude a given oasis. 

• For the near future the correlations involving S^^ and B{Bs,d — ^ /^^/^ ) will be 
most interesting as the data on these three observables will be improved in the 
coming months sharping the outcome of our analysis, possibly ruling out some 
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oases and scenarios of the couplings. The plots in Figs. [3| [19] and [7] will be helpful 
in monitoring these developments. 

Of particular interest will be the study of the effects of right-handed currents. 
Here the recent constraints on the Wilson coefficients of primed operators from 
b — )■ sfi'^fi~ transitions had already impact on our results. In the future an 
important role in testing RH currents will be played by 6 — >■ sui/ transitions. The 



plots in Fig. 21, 22 and 23 exhibit the power of these decays in this context. 



While in Bg^d decays for Mz' > 5TeV Z' effects are predicted to be small, in 
— )■ TT+i/z/ and Kl — )■ n^uu decays they can be important in LHS and RHS 
scenarios even at M^' = 10 TeV and slightly larger scales. This is seen in Fig. [10 



On the other hand as seen in Fig. 16 in LRS scenarios for Mz' > 5 TeV it will be 



difficult to identify NP in these decays. 

We have demonstrated that the imposition of U{2)^ symmetry on the Z' couplings 
has a profound impact on the correlation of S^^ with — )• /i^/x^ sharping the 
predictions significantly. We show this in Fig. 



25 



Our analysis of flavour-violating Z-couplings shows that in the case of Bd and K 
system they could constitute and important portal to NP with only small effects 
still allowed in the 5, system. The plots in Figs. |27l[28l[29l[30l|3lJ[32l[33]and[34 
illustrate these findings. 



We have reemphasized, following 44 46 , that large NP effects in rare K decays 



could be softened by the correlations with e'/e if the hadronic matrix elements 
relevant for this ratio were better known. 



We are aware of the fact that some of the correlations presented by us would be washed 
out if we included all existing uncertainties. Yet, our simplified numerical analysis 
had as the main goal to illustrate how the decrease of theoretical, parametric and 
experimental uncertainties in the coming years might allow to exhibit certain features 
of NP, even if deviations from the SM will be only moderate. In this manner we 
have uncovered a world of correlations present in NP scenario, where new effects are 
dominated by flavour- violating couplings of a heavy neutral gauge bosons and/or SM 
Z boson. In fact within the coming years the size of the assumed uncetainties in our 
analysis could likely become reality. 

We are looking forward to improved experimental data and improved lattice calcula- 
tions. The correlations identified in this paper will allow to monitor how the simple NP 
scenarios discussed by us face the future precision flavour data. 
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